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Reference input elements 
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I, 1 = ideal value 
Ye» Ye = system error 
Q, gq = indirectly controlled variable 
E/R 
C/R 


i 


actuating signal ratio = i1/(1 4-G H) 


control ratio = G/(1 +G H) 
B/E = 1100p ratio = GH 7 
B/R = primary feedback ratio = GH/(1 + G H) 
Where G is forward loop transfer function. | 


W, w = overall system transfer function or closed 
loop transfer function 


aWyi=2£W¥ 
RHP, LHP = Right or left half plane 
Large letters represent frequency domain (or transforms) 


and small letters represent time domain. 
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ABSTRACT 


This thesis contains a compilation and comparison of 
some of the techniques in use today for synthesizing and 
analyzing servomechanisms. The techniques are divided 
into the categories of Trial and Error Design Methods and 
Analytical Design Methods. Trial and Error Design Methods 
include such items as Root Locus techniques, Bode and 
Nichols plots, and other frequency plane plots. 
Considerable space is devoted to the Mitrovic method and 
to the Ross-Warren/Mariotti technique of compensatione 


Analytical Design Methods are generally considered 
to be those which use some definition of error as a 
performance index, where the objective is to minimize (or 
maximize) the performance index. Considerable space is 
allocated to those methods published by Newton. Some 
effort is also devoted to the problem of obtaining a 
proper statement of specifications. 


Thesis Supervisor: Robert Ke. Mueller 
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Chapter 1 
THE PROBLEM 


1.1 General. 

Feedback Control System synthesis means the 
determination of system and component specifications 
to meet the requirements of a specified job. The first 
phase of such a procedure involves the selection or 
design of a power element adequate to drive the load. This 
paper, however, is not concerned with the selection of 
components, and the power element is stipulated. 

The type of Feedback Control System considered is a 
servomechanism. By definition here, a servomechanism is a 
particular type of feedback control system in which the 
controlled variable is a mechanical position. The out- 
put is the mechanical position of one object relative to 
another. 

The words "synthesis" and "analysis" are frequently 
used interchangeably and loosely. Here, we will use 
J. Re Burnett's?” definition: 

(1) The synthesis problem. Given the input to a 
system and the required output, determine the transfer 
function of the system. 

(2) The analysis problem. Given the input to a 
system and the transfer function of the system, find the 


* numerical superscripts refer to Bibliography 





output of the system. 

The type of system under consideration is lumped 
parameter, finite, time-invariant, and linear. 

Stability is a problem to be reckoned with. We can 
say that stability is the primary consideration in all 
control systems. In that sense, it can be considered as a 
basic specification, always implied. However, in another 
sense, it is of only secondary consideration, because it is 
virtually always possible to render a system stable by some 
form of compensation. For linear systems, stability is a 
function of the system alone and is not dependent upon the 
input to the system. A system whose response will 
eventually become arbitrarily small, once the input is 
removed is "stable". Stated another way, a system is 
stable, if the impulse response of the system approaches a 
constant value and remains constant for large values of 
time after the impulse has occurred. If the closed-loop 
system function is stated as 

f(s) Pe 374+Gn., = Sa ee = NG) 
R(s) SGeba ae D(S) 
then D(s) is a Hurwitz polynomial, all of whose roots mst 


(m>n) 





have a negative, non-zero, real part, for stability. 

In the preparation of this paper we have researched 
a great many writings. We find all of them in agreement 
on one thing: In the design of Feedback Control Systems, 


arriving at an accurate and complete statement of the 





problem is the most difficult phase. The difficulty of 
this phase is aggravated by the fact that Feedback Control 
Systems have permeated all walks of our modern life. The 
science of feedback control knows no bounds as it cuts 
across all historic lines of endeavore The science has 
become so broad with so many people affected by it, that 

a statement of a particular design problem has little 
universal application. Each problem statement seems to be 
unique unto itself. 

The vastness of this new science is evidenced by the 
fact that many engineers no longer call themselves 
"Electrical" or "Mechanical" or "Aeronautical" Engineers, 
but rather "Control" Engineers. Many books are being 
written and published as "Control Engineering" books. 

There is a great deal of discussion encouraging the divorce 
of Feedback Control Systems from all other historical lines 
of engineering, to let it grow as an independent science. 
le2 Flow Diagram of Design 

A typical design flow diagram might be as follows: 


(See, for example, writings of R. B. Wilcox’). 






Problem 





Assigned | Statement; 





Note that "feedback" continually exists through all the 
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design phases. This is what makes control system design 
so fascinating, but at the same time so maddening. 

The problem statement may be given in any one of 
three forms: 

(1) a literal description of what the system is 
supposed to do, or 

(2) a numerical specification of the system per- 
formance, or 

(3) a graphical representation of the system. 
From the designer's viewpoint, numerical specification 
is preferred to eliminate ambiguity, but the price is 
loss of flexibility in meeting changing design considera- 
tions. If no graphical representation is provided the 
designer is immediately required to prepare a functional 
block diagram of some kind so that he may proceed with a 
component design analysis. Completion of this phase leads 
to the determination of suitable transfer functions for the 
components. Component transfer functions may be determined 
by experimental or analytical means. There are two basic 
experimental methods; determine the transfer function from: 

(1) the component frequency response 

(2) the component transient response. 
Analytical methods require the use of basic mathematical 
and physical principles (for example, Newton's Laws of 
Motion). 

In any event, a system block diagram may then be 
constructed, from which the "system" (or what we prefer to 


+s) 





here call "plant") transfer function is determined in 
either algebraic or graphical form. A pitfall comes 
to light here, because at this juncture, there is a 
tendency to forget the limitations imposed on the 
assumed mathematical model. These limitations might 
result from approximation, initial conditions, noise, 
drift or some other special input. 

System analysis might then proceed using time domain 
or frequency domain methods with or without the aid of 
analog or digital computers. It is this analysis which 
then forms a basis for any synthesis or modification 
technique. 

1.5 Problem Statement 

The type or manner of statement of the problem 
requirements often determines the feasibility of any 
particular approach. Although ideally, the requirements 
are documented in the form of complete numerical specifi- 
cations, more often than not, the only information is the 
general requirement that some given operation must be 
automatically controlled. The designer, therefore, 
must conduct a study to formalize his own detailed specifi- 
cations. 

1.4 Component Design and Analysis 

This paper is not concerned with this phase of 

design. It will here be assumed that the components of 


the plant are specified. We are concerned with the 
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component design of the compensator only in the sense 
of whether the component is physically realizable. 
1.5 Transfer Functions 

Since we here assume linear time-invariant systems, 
transfer functions may be manipulated following algebraic 
rules; isolation between component transfer functions 
is implied. It follows that the concept of superposition 
is valid. This means the output of the system may be 
determined as the result of all the various inputs and 
disturbances, each applied separately. It is also 
possible to distinguish between dynamic and steady-state 
performance in each of these casese The error may also 
be analyzed both dynamically and in steady state. 
However, a transfer function represents the assumption of 
some mathematical model, and as such, the limitations 
imposed by the assumptions must be kept constantly in 
mind.e In the mechanics of attempting system optimization, 
for example, the limitation of saturation is always 
present. As the signal approaches saturation, the assump-= 
tion of linearity becomes less and less valid. 
1.6 System Analysis and Synthesis 

The time and frequency analysis of the system 
characteristics is an analytical and/or graphical predic- 
tion that the system will perform as the dynamic specifica- 
tions require. If it does not, then synthesis is 
required to determine necessary modifications and 


revisions to fulfill the response and accuracy 
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specifications. Synthesis and Analysis then go hand in 
hands; the final analysis is the one which defines all the 
characteristics of the completed designe We generally 
refer to system modifications and revisions as system 
"compensation". Compensation will mean the inclusion of 
that network which is necessary in order to: (1) make 
the system stable, and (2) meet the specifications. 
A study of the literature reveals this phase to be the one 
where writers and designers reach a parting of the wayse 
They proceed down one of two broad paths. The best 
descriptive titles for these paths that we have found are 
those used by Newton: 
(1) Trial and Error Design Methods, and (2) Analytical 
Design Methods. Some of the more lucid writings in the 
current literature on these two philosophies are those of 
Thaler and Brown’ for Trial and Error Design Methods and 
those of Newton* for Analytical Design Methods. 
1.7 Objective of This Paper 

In this paper, we attempt to bring under one cover 
a synopsis of the current methods suggested or in use for 
the design of the servomechanism. Our objective is to 
provide insight into that phase of design surrounding the 
problem statement where the designer makes his decision as 
to what approach to attempt first. To that end an under- 
standing of the capabilities and limitations of the 


various techniques is required. The methods considered, 





although not all inclusive, are those which, in the 
opinion of the authors, do have the greatest usefulness 
or offer the most potential. 

In Chapter 2 we explore more thoroughly the statement 
of specifications. Chapters 3 and 4 contain a summation 
of the basic procedures and mechanics employed in the above 
two design methods. In Chapters 5 and 6, we have listed 
some of the modifications and refinements that are found in 
the current literature. The compilation is not all 
inclusive but the techniques listed attempt in some way to 
make the problem of compensation easier within their 


specified limitations. 





CHAPTER 2 
SPECIFICATIONS 

2el General 

It may often happen that the specifications given to 
the servo designer are either incomplete, incompatible, or 
incomprehensible. Specifications fall into two main 
categories: (1) The specification of control-system 
dynamics and performance, and (2) the general specifications, 
such as power-supply variations and environmental conditions, 
which influence the dynamic characteristics and performance. 
We mentioned earlier that the designer usually mist conduct 
a study to formalize his own detailed specifications. 
This study may immediately reveal contradictory specifica- 
tions. It may indicate that relaxation of a particular 
requirement will greatly simplify the control system. It 
may indicate that one specification is so domineering that 
satisfaction of that one is tantamount to overall success. 

A servomechanism is expected to perform any or all of 
the following functions: 

(1) Bring about a change in the actual value of the 
output so that it conforms to a desired value at all 
times; 

(2) Minimize the effect of varying component per- 
formance on the output; 

(3) Minimize the effect of disturbances. 


The performance specifications determine the degree of 
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excellence with which the servo must carry out the above 
functions. This means that the performance specifications 
must be given in terms of the desired output for a given 
input. It is, therefore, necessary to explore the type of 
inputs that might arise. 
2ee Inputs 
The input signal may be one of the following types: 
(a) Aperiodic, noise free 
(b) Aperiodic, with noise 
(c) Periodic, noise free 
(d) Periodic, with noise 
(e) Stochastic, noise free 
(f) Stochastic, with noise 
(Noise is regarded as any input-signal variation that is 
not a measure of the information carried by the input.) 
Some typical specifications for the six types of 
input signals listed are as follows (see chapter 16 of 
reference 4); 
(a) Aperiodic, noise free 
(1) The system dynamic error shall neither exceed 
a specified maximum value, nor shall the 
steady state error exceed a specified maximum 
value. 
(2) The integral-square system error shall not 


exceed a specified value. 


TE 





(3) The system error, in addition, shall not exceed 


a specified amount in the presence of a 
specified disturbance that occurs at some 


specified point in the system. 


(4) See paragraph 4.2 for other possible error 


specifications. 


(b) Aperiodic, with noise 


(a) 


Given that the first component of the system 
error is that for zero noise step input; the 
second component is the value of the output 
from noise alone; then the square root of the 
sum of the squares of the two components shall 


not exceed a specified value. 


(Aperiodic signals commonly considered are steps, ramps, 


impulses, pulses, or an input expressed as a power series 


in time.) 


(c) Periodic, noise free (only fundamental frequency 


present) 


(1) 


(2) 


The frequency response shall be characterized 
by a specified peak magnitude ratio (output/ 
input) occurring at a specified frequency. 

The magnitude ratio (output/input) shall be 
within a band of some specified number of 
decibels over a specified frequency range, and 
phase shift (output/input) shall not exceed a 


specified amount over this same range. 


i2 





(d) Periodic, with noise 

(1) Error expressed in a similar manner to that of 

aperiodic, with noise, abovee 
(e) Stochastic, noise free (see paragraph 2.4 for 
definition. ) 

(1) System error shall not exceed a specified rms 
value when the input autocorrelation function 
has a given value. 

(f) Stochastic, with noise 

(1) System error shall not exceed a specified rms 
value when the input signal autocorrelation 
function, the input-noise autocorrelation 
function, and the signal-to-noise cross- 
correlation function are given. 

2eo Disturbances 

It is also necessary to specify performance in 
response to a given load and/or disturbance occurring at 
points different from the input. The load or disturbance 
can also be classed as aperiodic, periodic, or stochastic. 
Typical specifications take on the same form as those above 
for inputse A load specification, however, usually 
prescribes the amount of time allowed for the output to 
recover to within a specified deviation. 
204 Stochastic Signals 

A stochastic process is one in which there is an 
element of chance. Sometimes the input to a system is not 


completely predictable and cannot be described by a 
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mathematical function of time. A typical example of a 
stochastic process is a radar signal mixed with noisé. 
Since the value of a stochastic signal cannot be determined 
with certainty at a given instant of time, probability 
density functions and other statistical characterizations 
such as the average value, the rms value, and the correla= 
tion function are used to describe the signal (see Appendix 
B). However, it is necessary to think of a stochastic 


2 each generated 


signal as a member of a family of signals, 
by an identical process. Such a family of signals is called 
an ensemble and the statistical characterization (such as 
a correlation function) of the stochastic process is related 
to the ensemble rather than to a particular member of the 
ensemble. It, therefore, follows that the determination of 
the response of a system to a stochastic input does not 
yield a function of time, but rather a statistical 
characterication of the output signal ensemble (see Chapter 
4) 
209 Philosophy on Choice of Test Input 

It may be necessary for the designer to prescribe his 
own inputs when analyzing the effects of noise, or load 
disturbances, or the effects of environmental conditions 
such as temperature, humidity, corrosion, etc. Even though 
the time domain characteristics are frequently specified in 
terms of the response ratio of the output to a step-function 
input, many others may be specified or ifmplied. Further- 


more, they may be extremely complicated—requiring 
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graphical description either in the tim domain or as a 
power density spectrum. It might be well here to quote 
the feelings of some writers in the field. It is to be 
noted that they are not all in agreement. For example, 
the following quotation (page 308, reference 5) succinctly 
expresses one viewpoint: | 

"The characteristics of a particular servo should be 
determined by the actual input, the actual uncontrolled 
disturbances acting on the system and the actual output 
requirements. It is clearly not sufficient to assume the 
input to be a step in displacement or velocity, nor is it 
sufficient to require only that the transient response be 
well damped and that the velocity lag be small .... In 
general, the actual input and noise as well as the output 
requirements need a statistical description ocoe.e One can, 
of course, conceive of specialized servo problems, in which 
the input is a displacement step and in which the require= 
ments are based on the transient response; but they are the 
exception rather than the rule." (words of Ro Se Phillips). 

On the other hand, quoting from the same reference 
(page 18): "The performance of a servo can also be 
specified in terms of its response to a step function. 
The procedure of experimentally and theoretically studying 
a servo through its response to a step-function input is 
extremely useful and is widely used for a number of 
reasons. The experimental techniques used in such 


testing are simple and require a minimum of instrumentation. 
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The characteristics of any truly linear system are, of 
course, completely summarized by its response to a step- 
function input. That is, if the step-function response 
is known, the response to any other arbitrary input signal, 
can be determined. It would be expected, therefore, and it 
is true, that with proper interpretation the step-function 
response is a powerful and useful criterion of overall 
system quality." (words of I. A. Getting). 

But, progress and advances in technology must be 
recognized. Less than ten years later, Dr. T. GC. Fry, Bell 


© said: ".... in any actual guidance or 


Telephone Lab, Inco 
fire-control problem, we are not really concerned with the 
response of the system to some particular, ideally defined 
tactical inpute We are concerned with its response to the 
whole gamut of possible tactical situations, including all 
the possible variability in the enemy path and all the 
possible errors which may be produced by random perturba- 
tions in the input data or in the mechanism itself. 
Obviously, this adds elements of information theory, 
statistical theory or whatnot and greatly increases the 
level of (essentially mathematical) insight required for 
effective work." 

The writers seem to be hitting at the same old saw of 
the tug-of-war between abstract mathematics and 
engineering approximations. To easily and practically 
accomplish a design by hand, on paper, requires simple 


inputs. The use of the step input, for example, does take 


6 





- aaa 


into consideration the the whole frequency range. But, 
as control systems become more and more sophisticated, it 
becomes necessary to consider more sophisticated inputs. 
This fortunately becomes feasible in the modern day with 
the use of computing machinery. 
2.6 Static Characteristics 

The static performance specifications describe the 
steady state value of the system output. Although the 
statement of them is a simple matter, their importance 
lies in the fact that they sometimes immediately establish 
the type° system which must be used. For example, if a 
system is to have no steady state error in response to a 
step input, it must be at least Type 1 (meaning one pure 
integration must exist in the open-loop transfer function); 
if no steady state error is to exist in response to a ramp 
input, the system must be at least Type 2. (Two pure 
integrations must exist in the open-loop transfer functions.) 

Normally an acceptable following error must be stated 
for the response to a ramp input. Some maximum allowable 
output is usually stated for the response to a disturbing 
input. 
2o7 Dynamic Characteristics 

The desired dynamic characteristics may be specified 
in terms of transient or frequency response. (If defined 
in both domains, their compatibility must be determined.) 
The three main themes of engineering design might be 
tabulated as follows: 


a) 
“Ns 





TABLE I 


SPEED OF RESPONSE | RELATIVE STABILITY 


Transient Frequency 
Domain Domain 












Transient | Frequency 
Domain Domain 


















a) Peak value — 










a) time a) bandwidth a)% of first 
constant overshoot of closed | 
loop output- 
input 
amplitude 





ratio 






























b) rise | b) cutoff b) number of gain and 
time frequency overshoots . phase margin 
of open=- 











c) settling] c) frequency of loop 
time peak frequency 
overshoot characteris- 





tic 





co 
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Allowable error in terms 






of % or per unit of controlled 







variable, stated as: 






a) Maximum 


b) Average 






c) rms value 
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2o8 Required Specifications 


The three main themes of design specifications are 
essentially those of (1) speed of response, (2) stability, 
and (3) accuracy. The Trial and Error Design Methods 
specification requirements are: 

(1) Input Signal 

(2) Desired Output 

(3) Disturbances and special inputs 

(4) Allowable error 

(5) Plant elements (fixed) 

(6) Relative stability 
The Analytical Design Methods specification requirements 
are $ 

(1) Input signal 

(2) Desired Output 

(3) Disturbances and special inputs 

(4) Performance index and required value of same 

(5) Degree of Freedom allowed in compensation 

Notice that the specifications for the two methods 
have two basic differences: 

(a) The Analytical Design Method calls for a 
performance index vice allowable error and relative 
stability specifications. A performance index is simply 
a Single number which is used as an indirect measure of 
system performance. Its use is an attempt to replace the 
functional description of the performance of a system 


through its response parameters (such as peak overshoot, 
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rise time, etc.) with a numerical description that rates 
the system performance with a single numbere 

(b) The Analytical Design Method calls for the degree 
of freedom allowed in compensation vice relative stability 
and plant elements specifications. This specification is 
not strictly necessary, but practically it is. With 
constraints imposed by degrees of freedom, it is possible 
to categorize large portions of the computational 
mathematics for all time, and make use of only the results 
of this categorization in a particular probleme More will 


be said about these differences in the next two chapters. 
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CHAPTER 3 
BASIC TRIAL AND ERROR DESIGN METHODS 


Oel General 

Almost all Trial and Error Design Methods depend 
upon the ability to express the output of a component or a 
servomechanism with respect to an input in terms of 
differential equations. These equations are almost 
universally transformed into algebraic equations by the 
Laplace transform. After this transformation is made, 
components can be collected into one overall mathematical 
model by well known, easily applied methods, and analysis 
and synthesis of the system is made directly in terms of 
these equations. 

Most methods utilize and depend upon the open-loop 
transfer function or equation, the closed-loop transfer 
function, or a combination of the two functions. Test 
inputs with simple Laplace transforms are applied to these 
equations, the output is compared to the input, and various 
parameters are measured to determine the acceptability of 
the mathematical model of the servomechanism. Then the 
necessary hardware is determined to match the model, or 
the model is changed and another comparison or analysis is 
made o 

Specifications and design parameters which must be 
satisfied by servomechanism adjustment are expressed in 


terms of either transient response to a nonperiodic test 
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input or frequency response to a periodic test input 
(almost universally a sinusoidal function). 

Several design methods work with the transfer 
functions directly in the frequency domain, taking 
advantage of the fact that the analytical or graphical 
results will be directly related to the frequency response 
parameters such as bandwidth or amplitude of response at 
resonant frequency. 

In the analysis of a servomechanism, most methods 
treat the servomechanism as a modification of a system 
with a second order differential equation. The output of 
any second order system is exactly known for any of the 
useful test inputs and the family of curves is not 
difficult to reproduce. This is done to simplify the 
description of the output with a given input to one that 
can be easily formulated. It is found in practice that 
the above treatment is almost always a useful approximation 
if it is but remembered that it is an approximation. 

Soc The Mathematical Model and Graphs of Response 

An early task of the designer is to assemble his 
likeliest components, or the given components, into a 
representative mathematical model. The whole process of 
design depends upon the validity of the model, being only 
as accurate as the model which is used to describe the 
assembly of hardware. The transfer functions of the 
components can be gotten from experimental tests made upon 


them or from the differential equation of the physics of 
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the component. Even when the form of the differential 
equation is known, experimentation is often required to 
obtain the values of the parameters. In this paper, 
analysis is described as the determination of the output/ 
input characteristics from the transfer function, but in the 
determination of the equations for a component all of the 
analysis theory and methods can be used in reverse. 

There are graphical displays which are used to portray 
either the mathematical model or the response of the system 
to a given input. Some of the more useful ones in analysis 
and synthesis are: 

(a) Frequency Response Graphs-=- amplitude of the output/- 

input ratio and phase shift vs. frequency. 

(b) Bode Diagram or Attenuation Diagram--log amplitude 
of the open loop output/input ratio and phase shift 
vse log frequency. 

(c) Nyquist Diagram- a polar plot of amplitude of the 
open loop output/input ratio and phase shift with 
frequency as a parameter. 

(d) Inverse Complex Plane Diagram--a polar plot of the 
inverse of the amplitude of the open loop output/ 
input ratio and phase shift with frequency as a 
parameter. 

(e) Nichols Chart= a log amplitude of the open loop 
output/input ratio vse phase angle with frequency 
as a parameter. Log amplitude of the closed loop 


(unity feedback) vs. phase shift is overlaid on 


oO 


the chart. 

(f) Root Locus Diagram-~ the complex plane upon which 
the poles and zeros of the mathematical model are 
plotted. 

O05 Specifications 
The operation of a servomechanism can be described 

by its specifications. These specifications are defined 

in such a way as to describe the speed of response, the 
stability, and the accuracy of the servomechanism. The 
specifications used in Trial and Error Design Methods all 
stem from the ratio of the output of the servomechanism 
to the input for a given type of test input. 

There are three major types of specifications. 

There are those specifications which are determined from 

the open-loop frequency response; others which are 

determined from the closed-loop frequency response ; and 

a third type which are determined from the transient 

response to a specified input, usually a unit step 

function. 

Theoretically the closed-loop response to a step 
input and the closed-loop frequency response can be shown 
to be equivalent * eo in practice, the correlation between 
the two is usually quite remote, and the conversion from 
one response to the other involves graphical methods of 
integration with many repeated calculations. There is 


direct correspondence between the open-loop and the 
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closed-loop frequency response and this relationship is 
not hard to determine. 

There are a great many different measurements used 
today in servomechanism design to specify the character 
of the servo. Most of these measurements are inter-~- 
related, and many of them are either synonymous or at 
least they describe the same type of output. 

The following list of specifications is intended to 
be mrely a sample of those more frequently ageds “eee 
Many of them require amplifying modifiers not to be 
ambiguous. 

(a) Open-Loop Frequency Response Specifications 

(1) Phase Margin is 180° minus the angular 
difference in phase between the output and 
the input at the highest frequency where the 
output is of the same amplitude as the input. 
It is a measure of the relative stability of 
the system, being unstable at non-positive 
Values. 

(2) Gain Margin is the ratio of the amplitude 
of the input to the amplitude of the output 
at that frequency where the phase shift 
between the input and the output is 180°, 

It is usually expressed in decibels and is 
a measure of relative stability, being 


unstable for non-positive values of decibels. 
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‘S52 (1) Phase Margin 


(2) Gain Margin 


(bd) Closed~Loop Frequency Response Specifications 


(3) 


(4) 


(5) 


(6) 


Peak Amplitude Ratio (Mp) is the maximum 

ratio of output to input. It is a measure 

of the relative stability of the system. 
Frequency of Peak Amplitude Ratio is a 

measure of the speed of response of the system. 
Bandwidth is one of the specification parameters 
which has no standardized definition. It is a 
measure of the frequency range at which the 
amplitude falls within specified limits. A 
popular limit is that the amplitude ratio be 
between + 3db and -3db. It is a measure of 
speed of response. 

Cutoff Frequency is much like bandwidth. It 

is the upper frequency at which the amplitude 
ratio reaches some specified value. Common 


ones are 0 db, ~ddb, =-6db, or -20db. 
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(3) Peak Amplitude Ratio 


(4) Frequency of Peak 
Amplitude Ratio 


(5) Bandwidth 


zi (6) mS (6) Cutoff Frequency 


(7&8) Damping Ratio (fj) ana Undamped Resonant 


Frequency (Wr) are specification factors 
in that they completely specify the 
response of a linear second order servo- 
mechanism. All of the other specifications 
listed are fixed by Ow and (Wn) for a 
second order system. The vast majority Om 
servomechanisms have a dominant factor 
which is second order in character. 
Consequently (7?) and (Wn) specify the 
type of response of a higher order system 


to a more or less close degree, 


(c) Transient Response Specifications to a Unit Step 


(9) Peak Overshoot (Mp, ) is the maximum amplitude 
of the first overshoot measured from the 
final steady state output. It is a measure of 
relative stability; the more the overshoot, the 
less stable the system. A system with no over- 
shoot is said to be "overdamped" while one with 
an overshoot is said to be "underdamped", 


(10) Time of Peak Overshoot (t,) is the time from 


ret 


(11) 


(12) 


(13) 


(14) 


step input until the time of the maximum 
amplitude of the first overshoot. Itis a 
measure of the speed of response of the 
systeme 

Rise Time is the time for the output of the 
system to cancel a certain portion of the 
errore A commonly accepted Rise Time is the 
time for the system to move from 10% of the 
final value to 90% of the final value. It 
is a measure of the speed of response. 
Characteristic time (€. or ¥Wnr) is the 
logarithmic decrement of the dominant 
portion of the response. It is a measure 

of speed of response. 

settling time is that time required for the 
servomechanism to reduce the error below and 
remain less than a certain percentage of the 
input step. Commonly accepted percentages 
are 5, 2, and 1 which makethe settling time 
become approximately 3, 4, and 5 times the 
characteristic time. It is a measure of the 
speed of response of the system. 

Number of Overshoots is the number of over- 
shoots of the output during the settling 
time. It is a measure of the relative 


stability of the system. 
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(15) The Error Coefficient is a measure of the 
steady state error of the system. It is 
dependent upon the type of the system and 
the gain of the system. 
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Oe4 Analysis 


A very important part of the synthesis of a servo- 
mechanism is the analysis of the system under test to see 
if it does or does not meet the required specifications. 
As we have seen, these specifications are invariably an 
expression of the comparison of one form of outputeto- 
input relationship or another. The more quickly one is 
able to determine this relationship, the shorter the job 
of synthesis will be. Many of the analysis methods 
available are in reality shorthand techniques at 


approximating the output-to-input relationships. 
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Remember that linear servomechanisms have output-to-= 
input relationships which are expressible as linear 
differential equations. The type of input to which these 
systems are subjected are the boundary conditions for the 
differential equations. Since it has been found that the 
equations can be handled or solved more readily in the 
Laplace transformed condition, mainly because the operators 
can be manipulated algebraically, the system equations are 
normally in this form. This is known as the frequency 
domain, and it is the one in which most synthesis is per- 
formed. 

The frequency response of a system is easily obtained 
by simply solving the closed-loop transfer function for 
jW. If the open-loop transfer function is known and 
available, rather than the closed-loop transfer function, 
one can obtain the frequency response directly without 
first solving for the closed-loop transfer function. This 
is done by plotting the frequency response of the open-loop 
system and then transforming the coordinates. If there 
happens to be a feedback function other than unity, the 
transformation requires one more step, but it is still 
useful. 

Unfortunately many of the specifications are pre- 
scribed in the time domain using the transient response to 
an aperiodic input. This situation has come abqut largely 
because it is easier for the specifications writer, and 


indeed for anyone, to visualize the effect upon, and the 
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action of, the servomechanism in that plane. So there is 
a requirement to know the response of the system to a 
transient input. There are several avenues of approach 
available to determine the transient response. 

The most obvious approach, but unfortunately the one 
which usually entails the most labor, is actually to solve 
the differential equation with the prescribed boundary 
conditions. Normally this amounts to taking the Inverse 
Laplace transform of the closed-loop transfer function 
multiplied by the Laplace transform of the input. A 
modification of this technique is to get the approximate 
output response by solving only for the dominant features 
of response, ignoring the less important features. Note 
that some prior experience is helpful in determining what 
is a less important feature. 

A second avenue of approach is to determine first the 
frequency response and then plot the time response from the 
frequency response. It has been shown that the two 
responses are uniquely equivalent. It is extremely un- 
fortunate that the equivalency is too obscure to be seen 
immediately by the average eye--or for that matter by many 
a trained eye. The transformation from the frequency 
response to the transient response to an impulse or step 
involves repeated graphical integrations, and without a 
digital computer can be a tedious operation. 


A third, and the most frequently used approach is to 


one 


assume that the system is like a second order system and 
that the relationship among its transfer function, its 
frequency response, and its transient response are closely 
related to a second order system. The great majority of 
servomechanisms are dominated by one pair of complex 
conjugate roots possibly with a single dipole near the 
origin. If the system is assumed to act like a second order 
one with only the complex conjugate roots, full advantage 
can be taken of the known relationships of second order 
systems with respect to transfer functions, frequency 
response characteristic, and transient response characteris= 
tics. The results can be appropriately modified if there 

is a dipole near the origin. 

The frequency response can be determined from the 
transient response. This is of value when a physical 
component is available for transient response testing, and 
the transfer function is unknown. 

In summary, it is seen that the loop transfer function, 
the frequency response, and the transient response of a 
servomechanism each uniquely defines its characteristics. 
Analysis is then the method of jumping from one of these 
four descriptions of the system to another. 

Sed Synthesis 

The determination of a transfer function or a set of 
components which will comply to given specifications is 
synthesis. Often our choice of components is either 
limited or partially specified. This restricts our ability 
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to determine the transfer function or the remaining 
components. 

In synthesis as well as in analysis, heavy reliance 
on second order (or possibly 3rd order) approximations is 
required in many of the present mthods, expecially if the 
specifications are given in terms of the transient response. 

An analysis of the specifications will quickly show 
what type of closed-loop transfer function is required in 
order to remain within the specifications. The real trick 
of synthesis is to translate this closed-loop transfer 
function, either real or implied, into an open-loop function 
which can then be stated in terms of components. 

This transformation is not a unique one and, therefore, 
there is an infinity of combinations of components which 
will fall within the specifications. Unfortunately there 
is a much greater infinity of combinations which will not 
fall within the specifications. 

The designer is allowed a considerable amount of 
leeway in the selection of a transfer function or of the 
components he must use. This leeway is allowed partly 
because the specifications usually require the system to 
fall within certain fairly broad limits and partly because 
the change in specifications due to a change in servo- 
mechanism parameters is usually not a sharp change. Since, 
to improve upon one specification, the servomechanism must 
often be alilowed to relax another specification, and since 


it is unlikely that the specifications can be described as 
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optimum to begin with, the design of servomechanisms is as 
much an art as a science. 

The most usual case is the one in which the whole 
system, except for a compensator or two, is already chosen, 
and it is only a matter of deciding what the compensator will 
bee Even the selection of the place to put the compensation 
may be restricted due to the limited physical accessability 
of the signal flow. 

506 Compensation 

There are two major methods of compensating a servo- 
mechanism. One is to place the compensator in the forward 
path (called cascade compensation), and the other is to 
place the compensator in the feedback path (called feedback 
compensation). 

Cascade compensation is the most widely used of the two 
because it generally requires simpler and less expensive 
elements. It is easier to synthesize a system using 
cascade compensation because the relationship between opene 
loop and closed-loop frequency responses and transfer 
functions is more direct. Drift of the parameters of the 
active elements in the forward path disturbs the effect of 
cascade Son wera ion more seriously than it would disturb 
feedback compensation. Cascade compensation elements are 
normally placed in the forward path at signal power levels 
rather than at output power levels. This allows the use of 


smaller network elements since little power must be 
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dissipated. The two types of cascade compensation that are 
generally used are lead compensators and lag compensators. 
Lead compensators act as a derivative path in parallel 
with a direct path. The derivative path tends to increase 
the output to input ratio at high frequencies. At these 
same frequencies it decreases the amount of phase shift 
that would occur in the uncompensated system. The impor- 
tance of the lead compensator is that it tends to cause the 
180° phase shift to occur at a higher frequency providing 
a larger degree of phase margin than would otherwise be 
available. The main objection to lead compensation is 
that it causes a great attenuation of the forward signal, 
requiring a substantial increase in required gain in the 
power element in order to maintain the same steady state 
accuracy as before compensation. Also lead compensators 
tend to increase the bandwidth of the servomechanism 
making it more susceptible to high frequency noise. The 


transfer function of the compensator is of the form: 
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Lag compensators act as an integrating path in par- 
allel with a direct path. The integral path tends to 
increase the output to input ratio at low frequencies 
while the phase shift is increased at the same frequencies. 
The lag compensator is used to increase the steady state 


accuracy of the uncompensated system. The time constant 
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of the pole of the lag compensator must be made large 
enough that the accompanying phase shift does not impair 
the stability of the system, and yet it cannot be too 
large or the transient error will tail-off too slowly. 
The transfer function of the compensator is of the form: 


es + | 
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Lead-lag compensators complement the virtues of the 
lead compensator with those of the lag compensator. The 
relative stability is increased with the lead compensator 
while the steady state accuracy is maintained, or at least 
the attenuation of the lead compensator is offset, with the 
lag compensator. 

Combinations of lead and lag compensators can be used 
to gain further benefits from these compensators. It must 
be remembered that isolation must exist between adjacent 
compensators for the transfer functions to be correct when 
multiplied. This isolation can be achieved with a buffer 
amplifier or cathode follower. 

Feedback compensation, while more complex and 
pPapenaanye in components, allows a greater flexability than 
cascade compensation. At reasonably 1 eae gains, it tends 
to nullify the effect of the forward components around which 
it is placed, making the whole appear like the inverse of 
the feedback function. Thus it is an ideal method of 
replacing an undesirable function with a desirable one. 

A forward component with a shifting or uncertain gain or 


parameters in the transfer function can be made quite 
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rigid by the proper use of feedback compensation. 

Tachometer feedback compensation is a very common 
method of dampening a servomechanism to make it more stable 
or to reduce the transient overshoot. The steady state 
accuracy suffers with tachometer feedback, but this may be 
corrected by placing a filter in cascade with the tachom- 
eter in the feedback loop. 
oof Summary 

The Trial and Error Design Method might be summarized 
as follows: 

(a) The given specifications are: input signal, 
desired output, disturbances, allowable error, plant 
elements, and degree of stability required. 

(b) On the basis of experience or preliminary 
approximations, select a form of compensation. 

(c) Establish the parameters for the compensation 
exclusive of system gain. 

(d) Adjust the system gain in accordance with the 
stability requirements. 

(e) Analyze system to see if the error is satisfactory. 

(f) If the error exceeds allowable limits, repeat the 
process, using different compensations continue until error 
specifications are mt. 

Although the theoretical design objective is the attainment 
of the "best possible" servo, a more mundane design objec- 


tive is the attainment of an "adequate" servoo 
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CHAPTER 4 
BASIC ANALYTICAL DESIGN METHODS 
4.1 General 

The authors received almost all their training and 
schooling in servos using Trial and Error Design Methods. 
The particular problem that motivated them for this type 
of thesis was the fact that they were unable to recognize 
an inconsistent set of specifications. Background training 
was more than ample to select some method and to forge 
ahead to a design; if it failed to satisfy the specifica- 
tion, to start over; to keep starting over until a solu- 
tion was found or patience was exhausted, never really 
knowing whether a solution existed, or in the case where 
a solution was found, whether it was the best solution. 

Analytical Design Methods purport to solve this 
dilemma. And this is true in the sense that if it is 
assumed that a Derromenrace index is able to incorporate 
all the specifications into it, the method, being a pure 
mathematical minimizing or maximizing process, is able to 
immediately reveal, solely by the mathematics, whether a 
stated performance index can be achieved. 

Recall that the three main themes in specifications 
are speed of response, stability and accuracy; it seems 
reasonable to believe that some performance index that 
includes error, time and an appropriate weighting function 
for the error, can represent the specifications, or at 


least most of them. 
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In selecting a performance index, we need an answer 
to the question of what kind of output is desirable for 
the servo. If there were no uncontrolled disturbances 
(or noise) our goal would be to make the output follow the 
input perfectly. But, generally, in the presence of 
disturbances, if the servo follows the input perfectly, it 
will also do a good job of following the noise. To 
establish a figure of merit or performance index, then, 
requires a compromise. The performance index must be 
practicable, not too difficult to apply, and of general 
applicability. It should also be a measure of the average 
behavior of the servo, rather than be affected by short- 
lived deviations from the mean, or shifts in the time axis. 
Quite a number of performance indices are in use or have 
been proposed, some of which are discussed below. 


4.2 Performance Indices for Transient Signals 


(a) Performance Indices which are not time weighted: 


a 
Qa) R= (el)dt (called Control Area) 
© 
(2) te (called Integral 
8 ={ j2d)dt Absolute Error = IAE) 


DO 
(2) ZA 4 (called Integral-Square 
oor { e*)dt Brror - ISE) 


Two of the above indices (P,; and Pz) can be used in purely 


analytical procedures whereas Po contains discontinuities. 
Pj was proposed by T. Mo Stout’, Py by Fickeison & Stout®, 
and Pz by Hall’ and Sartorius? ®. Pz is sometimes called 

the Hall-Sartorius criterion. It can be seen that all of 


these indices heavily weight the error at the beginning of 
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the transient. Even though they favor rapid speed of 
response, they also tend to cause adjustments that are, in 
general, less damped than is usually required. From this 
viewpoint an IAE index is better than an ISE index. IAE 
is, however, normally treated with the use of analog | 
computers. Note that none of these indices, though, place 
a penalty on errors occurring late in the transient, and 
thus, there is no good assurance that the servo does not 
have a long tail off. To provide for this, and at the same 
time, recognizing that no physical system can respond 
instantaneously, another set of performance indices are 
time weighted. 

(b) Performance Indices which are time weighted: 


ba 
(4) a = \ +t eltdt (called Time Weighted 
Control Area) 


ba 
co) eae iv t |e (+)| dt (called Integral-Time- 
a Multiplied Absolute 
Error=- ITAE) 


eS 
(6) irs = fx e dt (called pines AEM 
ISE 


a -_ 
P, was proposed by Nims” and Pe by Graham and Lathrop?*, 


Again, Pg and Pe can be treated in purely analytical 
procedures, whereas Pe is usually treated with an analog 
computer. P. can also be treated with a set of what are 
called standard forms of which the foremost proponents 
are Graham and Lathrop. It is discussed further in 
Chapter 6. The major objection to these indices is that 


they tend to tip the scales the other way: they put too 
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much weight on the latter portion of the transient. They 
might be more practical, if they ignored error after it 
reached a certain limit. Unfortunately such an expression 
would again not be analytical. However, Ps is quite 
popular in the literature for computer use. 


(c) Another index proposed by Nimst! iss 
‘aw 
(7) RP, = i elt) dt | [{ Secelde] 


However, we can find no other literature, for or against it. 
It can be used where the objective is to make P, and Py, 

some value other than zero. Pz (ISE) is by far the most 
popular index for purely analytical methods. One reason 
for this is that error becomes more undesirable as it 
increases in magnitude. As the below figure shows, 


weighting the square of the error rapidly penalizes large 


Orror o Weight 






| a ron. spare 
weighting 
Error 
However, probably the greatest reason for the wide usage 
of ISE (and Mean Square Error discussed in the next 
paragraph) is mathematical convenience. There is a highly 
polished body of mathematical knowledge that has been 


developed particularly around the idea of a mean square 


value. 
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4.5 Performance Indices for Stochastic Signals 

It was indicated earlier that usual specifications 
concern the response of the system to typical inputs, 
whereas systems are actually subjected to random inputs. 
This is one of the basic reasons for interest in statisti- 
cal criteria. By far the most commonly used and described 
index is the root-mean-square criterion, first proposed 


by Ae Co Hall? 


and developed by James, Nichols and 
Phillips”. It is characterized by its mean square value, 


thusly: 


. T 
(3) R= e(t) = Li ailing er(b)dt 


T9700 2T 
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Another criterion suggested by Oldenburg and Sartorius! 


is characterized by the following integral5 


(9) = |€4)| = ae fel) |dt 


The use of the root-mean-square criterion was ine- 
spired by the writings of Wiener!* and it has been highly 
developed by many writers, of whom two good ones are 
Newton and Laning?°, 

Recall that determination of the response of a system 
to a stochastic input does not yieid a function of time, 
but rather a statistical characterization of the output 
signal ensemble. Stochastic signals consist of two 
classes? 

(1) The process is stationary if the statistical 


behavior of the process that generates the ensemble is 


independent of time. 
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(2) The process is nonestationary if the statistical 
behavior varies with time. 

Many books are written on the fact that if a stochastic 
process is stationary, it is possible to use a single member 
of the ensemble of the process to determine the process 
statistics. (Various books also justify the assumption that 
stationary signals are valid for use in servo work. No one 
has successfully "cracked the barrier" of non-stationary 
signals yet, except for special cases.) Furthermore, the 
average value of a signal from a stationary process can be 
found either by taking an ensemble average at a particular 
time or by taking the tim average of a single member of the 
ensemble. Because of this (see Appendix B), the mean-square 
value of a stationary signal is precisely equal to its auto- 
correlation function, evaluated with the argument equal to 
zero. If Pg is then selected as the performance index, an 
incredibly powerful tool results. 

The rms value of a random quantity can often be 
calculated practically under conditions in which many of its 
other statistical properties cannot. For example, if a 
random quantity has a normal probability distribution, all 
statistical properties are determined from its rms and mean 
value. 

It should be stressed, though, that the rms value of the 
error does not characterize the error completely, and use of 
this criteria may result in overlooking some important 
aspects of the problem. “e@ (+) is independent of the 
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distribution of e(t) in the frequency spectrum. Since a 
system usually operates with other systems connected to it, 
the other systems will transmit e(t) in a manner that depends 
on their dominating frequency and on the spectrum of e(t). 
Gillel® indicates as an example the stabilization of an 
airplane carrying passengers. Resonant frequencies of the 
order of magnitude of the durationof a human pace should be 
avoided because such frequencies are the most adverse to 
comfort. It is just as important to consider the frequency 
spectrum of the error as its rms value. Nevertheless, there 
is widespread use of Pg in design analysis. This seems to 
be due less to its intrinsic worth as a criterion than to the 
convenience attached to its calculation. 
4.5 Configurations 
By establishing standard configurations, much of the 
mathematics can be worked out for one time and only the 
results are then necessary for design. The configurations 
used are’ 
(a) Fixed Configurations 

In this configuration, all the physical elements, 
including the compensation, are essentially fixed. Op- 
timization procedure then consists of the adjustment of a 
few free parameters (such as gain, or one time constant of a 
filter, or perhaps the ratio of a pair of time constants of 
a filter) so as to minimize the performance index. This 
technique ts frequently used in conjunction with the Trial 
and Error Methods where the form of the system has been 
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fixed by other considerations and only the best numerical 
values of the free parameters are sought. On the surface 
it appears easy to use since the only requirement is to 
differentiate the performance index with respect to the 
free’ parameters and to set the partial derivatives equal 
to zero. This technique is also easy to apply in connec-~ 
tion with the analog computer. 
(b) Free Configuration 

In this configuration, quite a high degree of 
mathematical sophistication is called for, because here 
the entire transfer function (or weighting function, if in 
the time domain) of the system is allowed to vary in 
minimizing the performance index. This can only be done by 
the use of the calculus of variations. This is one of 
Wiener’s big contributions, producing an implicit transfer 
function in the form of an integral equation known as the 
Wiener-Hopf equation. Solution of this equation is 
extremely difficult; Wiener conceived a process for solution 
called "spectrum factorization", Although this kind of 
procedure is very satisfying intellectually, it, in 
general, leads to a transfer function which is physically 
unrealizable 

(c) Semi-Free Configuration 

This can be thought of as the usual type design 
problem, and is the type that one first thinks of in the 
Trial and Error Design vernacular: All the plant elements 
are specified as fixed elements, but there is no constraint 
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on the type of compensation allowed. This technique again 
frequently leads to a compensation which is physically 
unrealizable, but it has a very real practical advantage: 
it provides a goal to aim for in compensating by Trial and 
Error Methods. In providing the goal, it frequently 
indicates the direction one must take in the selection of a 
compensator type. 
4.4 Rudiments of Analytical Design 

All of the methods use the overall or closed-loop 
transfer function. Adopt the following symbology (taken 


from Newton”) (Same as that on symbols page, but repeated 


here for convenience ) 3 Wl (Ss) : w(t) 
pia ae 
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Because the closed loop function is used, it sometimes is 
not too easy to manipulate back to the original block 
diagram. However, this disadvantage (if, in fact, it is one) 
is offset by the fact that these methods do not involve 
solving for the roots of the characteristic equation, unless 
they are desired in the final completed paper design. This 
is a very real advantage if the characteristic equation gets 
above an order of about five. 

Practically speaking, the application of these methods 
can lead to failure if the limitations of the mathematical 
model are not continuously kept in mind. They almost 
inevitably lead to non-linear operation because the per- 
formance index demands minimum error. The minimization of 
dynamic error calls for higher gain, leading to saturation; 
and it calls for cancellation, or néar cancellation of the 
plant elements by the compensator elements, leading to wide 
bandwidth. (For example, design of a second order system by 
these techniques will invariably lead to infinite gain which 
is not only physically impossible, but practically, undesir- 
able.) 

(a) Transient Input 

Let us first consider procedures for a transient 
input. Use of the ISE criterion has been highly developed 
by Newton & Gould”9*; use of ITAE criterion has been 
developed by Graham & Lathrop!” (discussed in Chapter 6). 
The ISE criterion is represented ass 

Ty = J ve dt 
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By Parseval's Theorem, the integral can be rewritten as; 


= a ("Yels) YeCs)ds 
—g re 


oy. 
a 27 
where Yg(s) is the Fourier transform of ye (t). If Yells) is 


a rational function, Iy can be written in the form: 

re r*c(s)eCs) 

1 = amj 2 BLS) DES) 
where C(s) and D(s) are polynomials in s. Definite integrals 
of this form have been evaluated and tabulated in terms of 
the coefficients of the polynomials (for example, see Appendix 
E of reference 2). Evaluation of the integral then gives: 
Iy= ty(?,B - +0 Pe) 

where P represents the free parameter. 


Iy can be minimized by solving the set of K equations of the 


kit ; 
orm Die r 
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Except for the simplest forms, solution of this set of 
equations can be difficult. Many times, the easiest approach 
is simply to plot ly versus P;, holding other parameters 
constant; the minimum value can be determined closely enough 
for engineering purposes. An important condition for use of 
the developed integral tables is that D(s) must have all its 
zeros (roots) in the left half plane. ‘This is also the 
condition for a stable system. Therefore, the parameters 

must not be allowed to take on such values as to cause D(s) 

to have zeros in the right hand plane. Otherwise the integral 


table is invalid and, besides, the system is unstable. 
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Note that another limitation exists above: Ye(s) must 
be a rational function. (It is not meant to cite this as a 
disadvantage compared to Trial and Error Methods, since 
those methods, in general, require rational functions also.) 
See Chapter 6, under Newton's Method, for a way around this 
limitations 

The above scheme can only be used with fixed configura- 
tion systems. For semi-free configurations with transient 
inputs, Newton uses translation functions. This is also 
postponed to Chapter 6. | 

(b) Stochastic Input 

Let v(t) and q(t) represent stochastic input and 

output respectively of a linear system whose weighting 


function is w(t). From the convolution integral we know that 


4(t) = 5 wht) u-(t-t,) dt 
q. (+ +O) = f wilt) ar (t+t-t2) dt, 


From the definition of autocorrelation functions (Appendix B), 


and from manipulating the above integrals, the following 


relation is obtained: Pe 
aq (0) = fu (i) dt, { wlts) Ry (t+t,-t,) dt, 


which expresses the autocorrelation function of the output 
in terms of the autocorrelation function of the input and the 
weighting function. 


In a similar manner? pe 
Pg (0) = J wlte) Gy (tt) dey 
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which gives cross-correlation between the input and the out- 
put in terms of autocorrelation of the input and the 
weighting function. If the correlation functions and the 
weighting functions are Fourier transformable, then 
$qq(9) = WES) WOES 
Pog (SJ = W (s) Fwy (s) 
The latter expression is frequently used to evaluate W(s), 
since it is the ratio of the cross-power-density spectrum to 
input~-power-density spectrum. 

From the block diagram, Up (t) es i(t)-@(t). The mean-= 
square-error is identically the value of the autocorrelation 
function of error, with? = 0: 

After suitable manipulation, the power density spectrum is 


obtained as: 
$,, (3) =B;.(5) - W4IS,,(- 8) W928, (5) + WOSWOE©) 


This equation is suitable for many variations. For example, 
if u(t) = Uz (+) + Up (whe re Vay represents the data or signal 
component and Un (4) the noise component, and we assume them 


to be uncorrelated, then: 


By (s) =(1-we) [fi- w-9 [BG +WE)WCs) Ban) 
(Obtained by letting }(s)=f, (sHB,,,(Sjand PS-S - (s)=B_,(S) 


in the previous equations.) 
Now, the problem, similar to that described for transient 


input is to find the area beneath the error-power-density 
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spectrum along the imaginary axis and then minimize this 
error. This can be done by evaluating and minimizing the 


following integral: 
i ae 
rps ( 
ue(t) = fy(= 7) Yulsids 
7 Rl 


Of course, the complexity of the mathematics involved has 
increased greatly over that for the fixed configuration 
transient signal. Although the error may again be minimized 
by setting the partial derivatives equal to zero, it is 
usually best to proceed with a series of plots of error 
versus the free parameters. 

(c) Stochastic Input, Free and Semi-free Configuration: 

Newton shows that the actual output can be eliminated 
from the error expression, and the error then is a function 


only of the input and the desired output as follows: 

__ bo dbo Deo 

Ye (+) e L.(d—2) wl) ary dt +{u tt.) dt, W (t2) YO, (, -t,) dt, 
es a Z deo = be) 


He then proceeds to prove that minimization of the error for 
a physically realizable weighting function requires sate= 


isfaction of this expression: 


[@mlta) fey (tthdta Pile) =o Fe $20 


This is the Wiener-Hopf integral equation in the time domain, 
of which W(t) is the implicit solution, whe re Uy (+) 
represents that weighting function that minimizes Ue (t) 0 


It is an understatement to say that solution of this 
equation is difficult. There are a few rare situations 
where its solution in the time domain is self-evident. 
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Newton (Chapter 5, reference 2) has a very good expose on 
the explicit solution to this equation in the frequency 
domaine Suffice to say that its solution consists of 
Wiener's "spectrum factorization" technique whereby a 
function is split in such a manner that one portion consists 
of that component whose poles lie in the LHP and the other 
portion, that which lies in RHP. By discarding the RHP 
portions, the explicit solution of W ais) gives a 
guaranteed stable transfer function. This solution in the 
frequency domain is much easier if the functions are Fourier 
transformable. Some authors have developed a solution 
whereby the whole derivation for optimum transfer function is 
conducted in the frequency domain. 

Since solution for a free configuration is somewhat 
academic, the equation has been modified to provide for the 


semi-free configuration. We quote that solution here: 
G.(-s) By: G) 


Wee, (<) = -{Gets) Ges Br@))+ 


[Ge(-3) Gs CS) [7 BF (2) 


Although this expression is quite ominous at first glance, 
a breakdown of the symbols helps: 
> 
Wem(s) = optimum cascade compensation for minimizing Us (4) 


Gr 


fixed or plant elements 


De 





eo 
— —_- 
= 


= any factor of Z\(s) which includes all the poles 


~~ 
ta 
~~” 
0 


and zeros of ZA (s) in the LHP. 
v(S)42 component of ral (s) which has all its poles in 
LHP such that Y (s) = B (sda has all its poles 
in the RHP. 
16) = ACS) [a*cs) = remaining factor of L\(s) 
which includes poles and zeros ot LA (s) in RHP. 
So, if the correlation functions are known, then solution of 
the above equation for W.,,(s) becomes a mathematical 
factoring problem, where stability is guaranteed. 
If the fixed elements are minimum phase functions the 


above expression reduces to; 


We Ss ect (3) [Wr ¢ (s) 
where W,,(s) = 45:9 /Euisf, [2.0 


For our purposes a minimum phase transfer function is one 





which has no zeros in the RHP. Now, W,,(s) is the solution 
for the overall system transfer function for minimum ye (t) 
where there are no fixed elements. Therefore, the overall 
system transfer function is independent of Gp(s) if Gp(s) 
is minimum phase. What this means is that when the desired 
ouptut of a system is equal to the input, the above equation 
will always call for a Woy(s) such that the overall system 
transfer function is 1 (thus giving a mean-square-error of 
zero) if the fixed elements are minimum phase. And this 
sounds reasonable, since the ultimate goal of any compensa- 
tion is te cancel all attenuation and phase shift caused by 
ake) 





the plant over the whole frequency spectrum, in order that 
the output will be an exact replica of the input. Even 
though this is the ultimate goal, we, of course recognize 
that physically, the goal is unattainable. Hence, this 
method in general leads to a compensating function which is 
physically unattainable. The above methods rather non- 
chalantly assumed that the correlation functions were 
available, whereas actually, the obtaining of a correlation 
function in analytical form is another barrier to be 
surmounted. 
4.4 Correlation Functions 
(a) Derivation from Theoretical Considerations: 
Appendix B lists some autocorrelation functions that 
can be derived. Because of this, they are popular 
for use in design. Furthermore, their use does 
not represent too great a departure from the world 
of practicality. 
(b) Derivation from Experiment: 
Any sophisticated or complicated design calls for 
the determination of the correlation function by 
experiment. Newton? presents an approximate 
numerical procedure for computing them from 
oscillograph traces. He also indicates an analog 
computer method whereby a stylus is manually 
moved over the oscillograph tracings. He states 
that currently the M.I.T. Servomechanisms 


Laboratory prefers a numerical procedure because 
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computers can be used. 

After determination by experiment, the 
functions are frequently not analytical and some 
approximation scheme is necessary if they are to 
be used in analytical equations. For example, 
reference 17 describes a scheme for approximating 
a correlation function by a series of damped 
cosine functions. If a function is approximated 
by any such scheme, generally the whole design 
problem is of such size that computers are 
necessary for the entire design problem. 

4.5 Summary 
The Analytical Design Method might be summarized as 
follows: 

(a) The given specifications ares: input signal, 
desired output, disturbances, performance index 
and required value of same, plant elements, and 
degree of freedom allowed in compensation. 

(b) Classify the problem according to free, semi- 
free, or fixed configuration. 

(c) For free or semi-free configuration, use an 
appropriately derived formla. For fixed 
configuration, express the performance index 
as a function of the free parameters; minimize 
the performance index by adjusting the param- 


eters. 
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(d) See if the compensation thus determined yields the 
required value of performance index. If it does 
not, the specifications cannot be met; if it does, 
practical realization may begin. 

The design objective is the attainment of the "best 

possible" servo to meet the specifications, using the chosen 


performance index as the criterion. 
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CHAPTER 5 
TRIAL AND ERROR DESIGN METHOD REFINEMENTS 


Sel Introduction 

The purpose of Chapter 5 is to aid a servomechanisms 
designer to find an appropriate method of approaching the 
solution of his particular design problem. To this end, a 
summarization of the techniques found profitable by the 
authors is made along with references to more complete study 
of the technique involved. 

The design problem, when determined by the designer, 
will contain three parts: (1) the specifications which 
must be met, (2) fixed elements to be included, and (3) 
any limitations or preferences in the type of compensators 
to be used. The open-loop transfer function or frequency 
response of the simplest possible servo containing the 
required elements should be analyzed immediately in order 
to decide the amount and type of compensation required. 
This is done by synthesizing a system which will meet the 
specifications and then, using the compensators and 
variable gains, constrain the actual servo model to be 
similar to the system first synthesized. The final, or 
possibly feedback, step is to analyze the final servo 
model to insure that it does indeed meet the specifica- 
tions. More care must be taken with approximations in this 
final analysis than in the preliminary analyses. 

Following a guide to design methods available, some 


of the less well publicized methods are discussed. Some 
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root locus techniques, Mitrovic's method, and some trans= 
formation methods are explained to the extent that they can 
be used with the help of this paper. Theory and develop= 
ment of the methods is left to the references listed, with 
only a hint as to the validity of the methods. A 

knowledge of the basic frequency domain and time domain 
methods is assumed. The best general review of the methods 


available is either "Ordnance Engineering Design Handbook" 


or "Handbook of Automation Computation and Control, 


Volume zien 


O02 Guide to Design 


The various methods to design a servomechanism can be 
listed in any number of ways: type transfer function used, 
domain used, the graphical plane used, analysis or 
synthesis, etc. None of these groupings is an aid to the 
designer when he is attempting to make his choice for a 
design problem at hand. The breakdown of methods, shown 
on the following pages, while redundant, does allow the 
designer to concentrate his attention to the methods found 
fruitful by the authors. It is by no means an exhaustive 


collection. 
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root locus 

Ag ae) 
1+G 

Mitrovic 

Bode diagram 


Nyquist 
criterion 


Bode diagram 


solve equation 
G( jw) 


root locus 
solve equation 
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1-¥ 
root locus 


second order 
approximation 


third order 
approximation 
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functio 
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16,25 
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20 
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pare Do4 


pare 50d 
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re sponse 


Limitations 


unity feedback 


unity feedback 


unity feedback 
unity feedback 





BC 4 


BD1 


BD2 
BDS 


BD4 


CB1 
CBe 


CBS 
CD1 
CD2 


DB1 
DB2 


DC1 
DC2 


DF1 


DF2 


DFS 


Me thod 


inverse Laplace 
transform 


Bode diagram 
root locus 


second order 
approximation 


solve equation 
W( jw) 


root locus 


Routh criteria 


Bode diagram 


second order 
approximation 


curve fitting 
Nixon's method 


Chestnut & Mayer 


Bode diagram 


second order 
approximation 


Floyd's method 
Guillemin's 
Stallard 
Nichol's chart 


Nyquist diagram 


Inverse polar plane 
plot 
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Reference 


Limitations 


20 


Xo) 


pare 5ed 
53516 


20 


16 
19, 20 


3 
Z, 45 16) lo 20 


20 
par. 5.7(&) 
pare 507(b) 


) 
5, 19 


pare 5.6(a) 
Pare 5.6(b) 


4 

Xe) must know 
feedback 

20 must know 
feedback 

20 must know 
feedback 





Me thod Reference Limitations 


FA] Bode diagram 3 

FC1 Chestnut & Mayer 1, 45 19 
diagrams 

FD1 Nichol's chart Xo) 

FD2 Nyquist diagram 20 

FDS Inverse polar plane 
plot 20 


(Note: paragraphs indicated above refer to this paper.) 


oed Root Locus Analysis 
The root locus method was first developed by Walter 


Evans!® in 1948. It has gained considerable popularity 
with usage. There have been many modifications of, and 
techniques developed about, the root locus. 

Basic uses of the root locus are the plotting of the 
open-loop transfer function poles and zeros on the complex 
plane and the plotting of the locus of the roots of the 
closed-loop transfer function on the complex plane using 
the poles and zeros of the open-loop transfer function. 
The plotting of the root locus is basic to the general 
method and is assumed to be known by the dedtener” °t??"° ‘ 

The limits of gain which will allow a system to 
remain stable, and the second order approximation of the 
relative stability for a given gain, can be determined 
from the root locus. 

Analysis of the root locus, from the position of the 
closed-loop poles and zeros, can be determined by the 
following methods which are listed approximately in order 
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of increasing difficulty, but also in order of increasing 
accuracy 
(a) Graphical Representation of Dominant Pole-Zero 


Locations Versus Selected Specifications. 
Elgerd and Stephens“1 have graphed the step input 


response for almost every combination of roots and zeros up 
to fourth order for selected distances from the origin and 
for selected conjugate pole distances. 

Abbott and Patton** have graphed the time for the out- 
put to first equal the input, the peak overshoot, the time 
of peak overshoot, and the settling time (all for a step 
input) for closed-loop systems containing a complex 
conjugate root pair (1) alone (second order) (2) with a 
zero (second order), or (3) a real root (third order). 

They further show that systems containing a complex 


conjugate root pair, a real root, and one or two zeros can 


be broken into partial fractions of the above three types 
and the transient responses then added. 

Both of the above techniques require the system to 
be normalized to make one of the roots or zeros unity. 
Elgerd and Stephens consider more types of systems and 
their article (Trans AIEE, Vol 78, Pt II, 1959) is more 
apt to be available, at this time. The method of Abbott 
and Patton requires essentially no interpolation due to 
their mapping of the responses on the complex plane. Both 
methods show, along with Chu’s (see below) that a real 
pole or zero loses most of its dominant influence when 


it is more than twice as far from the origin than one of 
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the complex conjugate roots. See Figure 1 at the 


end of this chapter. 


NS 
(b) Chu's Approximation Equations 





This method computes, either graphically or analyt- 
ically, the dominant terms of the transient response to a 
step input. From this he finds simple equations for the 
peak overshoot and the time of peak overshoot. Let 


- (s) = kA) _ Kk IT (st 2) 
: BO T+ 4) 


/ 
A (Ss) + Ang Bis.) ( 


Cane 
Wp ale 
Me = Z k Acs) é 


To =e LU a Bese) 


‘Tp 





_at 


Clty= kB ++ ZZ iain C cos, (Wot Ves.) 


(a) 





+ k Aw) pan 


Tt hecessary only 
Cc) iF T, S 3/Tp 
whe re; K = a factorable constant 
Ato = the product of the distances from all 


the zeros to the origin 


IS = the product of the distance from all the 


roots to the origin 


U/\ 
O 
a 


ow + {W, the location of one of the 
dominant complex roots 
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A is.) = the product of the distances from all 


of the zeros to Ss, 


; 
Bos.) = the product of the distances from all 
the other roots to Sse 


Yes) = Ang A (s,) -Ang Bis.) — Ang So 
Ang Acs.) = the sum of the angles from the real 
axis to the line from the zeros to 8s. 
Ang Bis,)- the sum of the angles from the real 
axis to the line from the other roots 


FO Soe 


Ane Se = the angle from the positive real axis 


Ac) 


to the line from the origin to S85. 


the product of the distances from all 


a4 


the zeros to GQ, 


the product of the distances from ali 


C8 


) 
S (a) 
the other roots to Gy . 


(c) Wheeler's Approximation Equations<* 


Wheeler's equations are similar to Chu's, except 


that he makes the substitution: 


kK Bis Ks.) = acre 
Q A (So) A (S ) 


le a = Clo.) = Final value 
(0) 
Thus, if one knows the final value of Ce) or Fe(sso he 


may plot the step response from the open-loop poles and 
the closed-loop roots. In any case, he can plot the error 


to a step. 
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where: Ven = the product of the distances for all the open- 
loop poles to S,. 
(d) Graphical Residue“? 
The graphical residue gives the true transient 
response to an input step. If there aie repeated roots, it 
will not work in the form given below but the above reference 


should be consulted; 


wet 
at) eps +). kA ¢ 








T real Bien) 
l A aire 
i me —— = aaa nm 
Bice WW 
(+54), imag — Ang A,r -Ang 5 — Ang a) 
24 


(e) Frequency Response from Root Locus 


The frequency response of the system can be gotten 
graphically for any givenW by measuring the distances from 


all the roots and zeros to jw on the imaginary axis. 
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M (w) = en =z Magnitude. of response. 


N (w) = & Ang Aljw) - 2 Ang By) = Boose shi t+ 


response 

oe4 Root Locus Synthesis 

Having found where the desirable locations for the 
dominant roots are, the designer must now be able to 
"maneuver" the system's roots into these favorable locations. 
This can be attempted by varying the parameters of the given 
system or by adding a compensator in cascade. 

(a) Effect of Varying One Parameter of Root Locus<* 

Dre Wheeler demonstrated that any one parameter of 
the open-loop equation may be factored out of part of the 
characteristic equation and then be used as the gain of the 
final root locus, thus showing the effect on the system of 
varying this parameter. After factoring out the variable 
parameter, a root locus plotter would be a definite aide 
This type of plotter is available in different degrees of 
sophistication“©, The real use of this technique would be to 
see the true effect of one variable parameter. 

(b) Placement of Lead or Lag Compensators“’ 

The Ross-Warren technique is a natural outgrowth of 
the root locus equations and of earlier efforts by Walters“8 
and Aseltine“%. The locus is found to a point where the 
dominant complex conjugate roots are desired and the phase 
angle is noted. The placement of the lead or lag elements 
is partially fixed by this phase angle because they must be 
placed so as to change that angle to -180°. The compensator 
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is further, and completely, fixed by stipulating that the 
error coefficient (Ky) be unchanged by the compensation. 
This stipulation is unnecessarily restrictive and is one of 


the major drawbacks of the method. 
Measure ° pa 
PC « Ox , Ge, Une , es Jie ) K, 


P = (2n-1)T- OG, where nm is such that |W| < 180° 
h = cot" (et p- ese f) 
Ze $= QOS sind /G, sin @ 
Dp = Wr sind/sin(o-¥) 


whe re 3 

Oc = the phase angle measured to the desired 
root location. 

Ge = the product of distances from the zeros 
to the desired root location divided by 
the product of the distances from the 
poles to the desired root location. 

Wn, = description of the desired root, in 


second order terminology. 


K = the gain of the uncompensated equation. 
Zz = zero location of the compensator 
i = pole location of the compensator 

s p = attenuation of the compensator. 


A modification of the above method allows the 


designer to approximate immediately, a lead filter with a 
Su 





zero, and to use other than the uncompensated K,. See 
Figure 2 at the end of this chapter with which one enters 
with the angle to be compensated and obtains the zero 


placement°° 


« The effect of a real pole can be approximated 
from the same graph by subtracting the angle found at the 
pole distance from the angle found at the zero distance. 

Mariotti has developed a similar me thoda©? for the 
determination of pole-zero location. The Ross-Warren 
technique, modified, requires successive approximations to 
place the compensator at the place with the best error 
coefficient, while Mariotti has graphed this placement so 
that it can be accomplished enamaeee ty. His work shows 
that the error coefficient increases as the compensator is 
moved further from the imaginary axis. 
5.5 The Method of Mitrovic?9°t 

(a) Theory 

The method of synthesis presented by Du$an 
Mitrovié is a different and potentially powerful approach to 
the solution of the synthesis problem. This method utilizes 
the characteristic equation in the unfactored polynomial form. 
fi. =eel. Ss .Gy eee eee AQ,StAo 

Separating the real and the imaginary parts of the 
characteristic equation and solving for a, and a1 produces 
two parametric equations in terms of the higher order 
coefficients and with the variables {f and U). Inspection 
of these equations reveals that each power of Ww contains a 
single coefficient of the characteristic equation and an 
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easily tabulated function of Ff ° 


a, = —w [a4 tas Awe + s+ Onth., (pw) 


a, = Ai. (P) uw + Qs, (3) w+ ao +n f(Dw” | 
, 2 OD 
#.CY) = “Jef Puy t+ Pea 9)| maa ier 
The basic consideration of the Mitrovit method is to 
let a> and a, become variable and to observe the (aj,;8) 
plane. 
(b) Properties 
The locus r( ¥) of the parametric equations is 
plotted for a specific f and is found to have some 
remarkable properties: 
(2) 8 faa Co) 
The servomechanism is stable if certain conditions 
are met. This is because the total phase rotation 
as (Ll) » OO is fixed for a stable system similar in 
nature to the Nyquist criteria. This curve 
immediately shows the range of a, and a, for which 
the system will be stable. 
(2) -(f,) 
The characteristic equation has no roots with f 
greater than Meat certain conditions are met. 
This curve immediately shows the range of a, and - 
a, for which the system will have agit less than 
4, o More important, it will show the exact 


aj] and a, necessary to make the system have roots 
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(3) 


(4) 


at any givenl) ., and vl Ke 

ri) 

The characteristic equation has all real roots if 
certain conditions are met. Here again, the range 
of a; and a, are shown. The real roots are 
quickly determinable. The effect of changing 

4); Or a) is immediately evident. If there is 

but one complex pole pair, it can be determined 
easily once the real roots are known. 

Third order equations 

If the system can be approximated by a third 
order equation (i.e. a dominant complex conjugate 
pole pair and a maximum of one dominant real root 
with no more than 3 dominant zeros) this method 
is considered the ultimate in flexibility by the 
authors. 

A third order characteristic equation can be 
normalized so that ag and ao are both one. Now 
all third order equations can be solved by one 
set of ry. This set is tabulated once and 

for ai1°°°+, ‘The 4, and a, required for all 
dominant root combinations is immediately 
discernable. Thus a very powerful tool is 
available to determine the parameters of any 
system with a closed-loop characteristic equation 


of third order or less. 
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(c) Plotting of aie) 


The equations developed in the theory can be used to 
plot rcs) with a table of presley). 

An easier method is available, though, using a 
table that has the ip (§) Waa constructed 
for the ei in question. Tables II, III, IV and V 
(at the end of this chapter) are constructed for 
values of if of 0, 0.5, 0.7 and 1.0. 


Qo = —W*[de Wl w) +e Vs (FW) 4 Ant 64 
tig On (Y,w) +43 %(L w)+ 10 t n V($w) 
YC YW) = Ge (F)wk! 


These tables are constructed for values of LA) 
from zero to one. To insure that the area of 
interest falls within the above UW), the characteris- 
tic equation must be normalized so that both ay, 
and @_ 1 are unity. This is accomplished easily 
by the substitution s = an_}7 Pp, and it insures 
(for equations with all negative real parts in the 
roots) that all roots lie between O and =-1. The 
transformed parameter is related to the actual 


parameter by: 


es, 
WwW. = —_ 


EL, 
ee = ~—- - | 
ot = “a * Cer (An-1)” 
(d) Determining Pret 
The [(o) curve is plotted and point (aj 5a,) is 


Jill 


noted. 


The system is stable if: 


(1) The point (a, 5a.) lies in the first 


quadrant, and 


(2) while (J varies between C and + 09 the 


(3) 


curve alternately cuts the lines a, and 


&a) with the provision that a. is cut first, 


re) 
and 

that the total number of points of inter- 
section is n, where n is the order of the 


characteristic equation. 


(e) Determining Y 


the T(¥,Jcurve is plotted and the point (a; 5a.) is 


noted. All roots have damping factor greater than Sie i 


(1) 


(2) 


(3) 


The point (a,, a,) lies in the first 
seiashaceeey, and | 
Sileucurvalmengiiwenaiel: Gaucuthenlinse 
&o and a, with the provision that a, is cut 
first, and 

that the total number of points of intersec- 
tion is m, where mis the largest integer 


which fulfills the inequality: 
2 © 
me [1-1] & g n(i+ a 


where n is the order of the characteristic 


=| 
equation and 2 = = — COS J é 


If the point (a,,8 5) lies on the f (Y,)curve , 


ve 





the damping factor of a complex conjugate 
root pair is {) and Wy, is the W of that 
point on the curve. 


(f) Determining the Real Roots 
The rq) curve is plotted and the point (a, 5a) is 


tae 


noted. 
Ai weal roots 12d -onnd f7(\) eienewde S Geek v7 
(when Ay = 1.0, W = O~) from the origin such 
that a line tangent to af will pass through (27585) « 
Conversely, at any place that a line from (a},a,) 


is tangent to the curve, there is a root on the 

negative real axis at Cees 
(g) Lead and Lag Compensation 

The use of the Mitrovic method is quite easy with 
lead or lag compensators. The characteristic equation 
of a lag compensated transfer function will probably 
have the lag pole in coefficients other than a, or a). 
Lt will be found that the effect of the lag pole on 
these coefficients is so small that it can be deleted. 

The characteristic equation of a lead compensated 
system has the same problem of having the compensator 
pole in the higher order coefficients. Unfortunately, 
it is a larger term and cannot be ignored. Instead it 
must be fixed, in order to get numerical vaiues for 
the coefficients. This would also have to be done to 
the zeros and the gain if they appeared in the higher 


order terms. Then any computations involving 
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variations of parameters which were fixed in the 

higher order terms would either be approximations or 

the curve would have to be continuously replotted. 
506 Frequency Response to Transient Response 

The several methods for obtaining the transient 
response to an impulse input from the frequency response 
curve make use of the following relationship; 

bo 
f= Ele ut RLFC) dw 
© 

where AA [F Gu) } is the real part of the frequency response. 

The real part of the frequency response can be computed 
directly from the closed-loop frequency response curve. If 
the frequency response were plotted in polar coordinates, 
the real part could be read directly. Chen and Shen? 
suggest a chart which has the closed-loop frequency 
response plotted in polar coordinates and overlaid upon it 
is the open-loop, unity feedback frequency response. 
Gould? recommends a chart, similar to a Nichols chart, on 
which the open-loop frequency response curve is plotted as 
log-gain vs. phase angle (See Figure 3 at the end of this 
chapter.) Overlaid upon the chart is the magnitude of the 
real part of the closed-loop frequency response for a 


unity feedback system. 
4,19 


(a) Floyd's Method ° 


The real part of the closed-loop frequency 
response curve is approximated by the sum of a series of 


trapezoids: 
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=p >i 
Oe 





A check on the it ate of the trapezoids is: 


£(0) -2> 1 Wi 


(b) Guillemin’s Method id 

| The straight line approximation of the real part 
of the closed-loop frequency response is differentiated 
twice, leaving a series of impulses of height Ci and at 
frequency WwW 

mM 
C+) 7 a Go cea ) a; Cos we | 
a a ie a : 


A check on the correctness of the impulses is: 


flo) = Pay 


(c) A Computer Haitiod 





Levadi has combined Floyd's method and Guillemin's 
method in a way that is easy to program on a computer. 
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The output from the computer is the impulse response, the 
step response, and the response to a ramp inpute 
oe/7 Transient Response to Frequency Response 

(a) Nixon's Method? 

A straight line approximation of the transient 
response to a step input is made and it is broken into the 
sum of a series of ramp functions. These functions could 
be expressed in the time domain as a magnitude multiplied 
by a unit ramp with a time delay. Instead, they are 
expressed directly in the frequency domain (by Laplace 
transform) and multiplied by se This is the transfer 
function because, before multiplying by s, this sum is 
the Laplace equation for the output to a unit step input. 
i the transient response to an impulse were used, the 
multiplication by s would have been unnecessary. To find 


the frequency response, _ jw for s3 


cio) = > (a E ie | i 


Since the equation is imaenaa at i) = = 0, Ne ee 


-~_ «= 


rule must be usede 
(b) Chestnut and Mayer's Method4 


The transient response to a step input is broken 
up into the sum of a series of step functions each occurring 
at an equal interval of time, but being of a variable 


rb 


height AC;: J 


fas. 





-_ 
= ‘ oe 
—_ _ ea: a= ee ye _ 
' 
_ 
— 
= 
a ee 
= — 
= = 
es «¢ 
= 


n -ty (fu) 
= () = 2 ACae 


The transient response to an impulse is broken up 
into the sum of a series of pulses of uniform width At 


but of variable height C4; 


n oe 
F(jw)= 2 Crate gw) 
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TABLE II 


TABLE FOR COMPUTING MITROVIC'S [(Q) CURVE 


wl A th Nl % [| | [% ho 


-l1| oI @ 
-1 
ol 
-1 
-1 
~1 
oJ] 
=} 
=~] 
oJ] 
-1 
-1 
-1 
=] 
=~] 
-1 
-1 
-1 
=] 
-1 
~l] 





O0000 0000000 00000000 
O0000000 00000 00000000 
OCO00000 0000 000000 
OCOo00000 0000000 0 
OCOo00000 00 0000 


WH) vs. p tor f=O 


YAY, ww) = Y,($) 0 
do = -W7Ja,wly w+ ds P.(4u) ++» + An Hho] 


G=Qv(¢ Ww) +4 As¥(¥ w)+ ot Gy YA (Y, w) 
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TABLE Il 


TABLE FOR COMPUTING MITROVIC'S [~ (0.5)cURVE 


w |e [vi)ve 4 bm [oe [Yel % [to |% | Ho 
0 O -1] 0 

















eO5 | .003]) -1].050] O | 0 

elO | .O1 -1/.100] O | 7001] O 

ol5 | .023}) -1]/-150] O | 7003] 001 

020 | 04 -1}-200} O | .008; .002 

e259 | .063| -1}].250] O | 3016] .004| O 

030 | 09 -1}.500} O |2027/] .OO08} O 

055 | o1235| -1}-550}] O | -043] .015} O 

o40 | 16 -1}.400]} 0 |3064| .O26] O O 

045 | .203| -1).450] O | 7091] .041] O 0 

050 | 625 -1/.500/] 0 | 2125 | .063} O 0 

09D | e503 | -1/.550] O | -166/ .O92| O O 

060 | 36 -1]/-600/ O |2216|{ 130} O O 

065 | e423) -1/-650|] 0 [7275] .179| O 0 

°70 | .49 -1|-700] O [7543] .240] O O 

of5 | 6963] -1]-750/ O | 3422; ~516} O O 

080 | .64 ~1 |.800 | O |; 3512] .410] O O 

085 | .723 |] -1].850) O | 2614] .522) O te) 

090 | 81 ~1/.900} 0 {3729} .656|] O re) 
0903 | -1/}.950 | 0 |}3857 |] -815] O O 
1.0 -1j1.00/ 0 | 1.0 160 | O 0 








Ww) VS, WY Se f=as 


(Sus) = eC $) Wr 
Ct, = -w*[|aa¥CF W) +A, C4, WH eee + Ay Yn (4, w) | 


OQ, = Qi (Putas XY wt ++ + On 4 (¥ wy) 
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TABLE IV 
TABLE FOR COMPUTING MITROVIC'S {(0.7) CURVE 


0.0 2-1; 0 0 






005 


32002 
27004 
2009 
7018 
3034 
2062 
7° LO8 
7180 
9295 
2463 
21135 
o LlL2rl.075 | 


Wvs.W for f=7 


VAY w) = & Cf) Ww 


Qo 


i] 


~w™ [G2¥(Y w) + og Y(fu)+ tanh (gua) 


G,= Us W, (fw)+s Y,C4 w)+ » ot An YW, (4 w) 
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TABLE V 


TABLE FOR COMPUTING MITROVIC's [(I.o) CURVE 





LU vs. Y for f=/.0 
4.9, w) = $$) 
Ty = -W* (O28 (Fw) +459, (Yue +a, ¥, Gui] 


Q,= C2 Bf w) +Gey,C¥ Ww) + ss +p W(L vw) 
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CHAPTER 6 
ANALYTICAL DESIGN METHOD REFINEMENTS 
601 Introduction | 

In this chapter we summarize selected methods from 
current literature that concern Analytical Design. Some of 
the methods do not strictly meet the definitions of the 
Analytical Design Method at the beginning of this paper, 
but they are included here, because the literature 
generally considers any technique which uses a performance 
index as its main criterion to be "analytical". 

6.2 Newton's Methods@ 

Chapter 4 was built almost exclusively around the 
writings of Newton, Gould and Kaiser“. However, they have 
done a great deal more than to develop some basic 
mathematics; they have expended much effort toward applying 
the basic mathematics. Their methods use ISE and MSE 
exclusively, since, as they state, these indices are the 
only ones of engineering usefulness that lead to reasonably 
straightforward mathematical analysis. 

(a) Translation Functions 
Newton proves, by means of variational calculus, 
that translation functions can be treated in a manner 
similar to correlation functions. Translation functions 
(see appendix B) are used with transient inputs. The 
form las indicated in Chapter 4, then, can be used by 
merely substituting the correct translation function in 


the place of the correlation function. Most translation 
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functions can be found by Parseval's Theorem. Use of step 
and ramp inputs, even though they have no Fourier transform, 
impose no special limitation, since it is possible to 
introduce a convergence factor to find a Fourier transform. 
(The Fourier transform found this way is the same as the 
Laplace transform; however, special care is required in 
taking the inverse transform). This scheme is also nuch 
easier to use if the fixed elements are minimum phase. 

What appears to be a particular advantage of the method is 
its indifference to the kind of transfer function that exists 
in the plant. For example, it is able to handle a 

=a Although this 
method normally employs the frequency domain for finding 


Gr (S)= ea! as readily as a Ge(s) = 


the best compensation, it frequently finds the minimum 
value of ISE by working exclusively in the time domain. 

(b) Saturation 

As indicated earlier, compensation by error 

minimization will almost invariably lead to saturation in 
a linear system, rendering the assumed mathematical model 
invalid. The reason for this is that the "called for" 
compensation transfer function is usually that which will 
cancel the effect of the plant transfer function. Since 
the plant normally consists of integrations and lags, this 
means that the compensation (or, at least, the equivalent 
cascade compensation) will call for several differentiations. 
Differentiations lead to large signal excursions in the 
system, with the result that saturation occurs somewhere. 


Differentiation also calls for more bandwidth, since the 
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Signals have high frequency components. Many physical 
systems have the characteristic of linearity for large 
amplitudes at low frequencies, but can only tolerate small 
amplitudes at higher frequencies if linearity is to be 
preserved. To make the Analytical Design Method practical, 
this reference introduces constraints into the equations to 
insure that saturation does not occur. 

(1) Transient Inputs - The most direct approach for 
avoiding saturation is to limit the peak value of 
that signal in the mathematical model that 
corresponds to the signal in the physical system 
that is likely to cause saturation. However, it 
is impractical to express the peak value of a 
Signal as a function of the free parameters for 
systems above second order. The reason is that 
finding the roots of the characteristic equation 
above that order requires the use of numerical 
values so specific values of the parameters are 
necessary. But the integral-square value of a 
signal can frequently be expressed in terms of the 
free parameters; since the integral-square value 
gives large weight to large values, then some 
degree of control of the peak value is obtained. 
One procedure is to express the integral-square 
error as a function of the free parameters; also 
express the integral-square value of the saturaration 


signal as a function of the free parameters. Then 
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solve the equations by adjustment of the parameters. 
For fixed configurations this is best done by trial 
and error. For free and semi-free configurations, 
this reference introduces a more refined technique, 
employing a method developed by LaGrange. By use of 
a synthetic function and a constant called the 
LaGrange multiplier, the minimization of error 
problem is converted from one with a constraint to 
one without a constraint. Either of these methods 
works well in a fixed configuration problem, such 
as a positional servo. In this type of servo, 
saturation most frequently results from demanding 
an acceleration of the output member or load, which 
the motor is simply not capable of providing. It 
is to be noted that each introduction of a constraint 
effectively eliminates one free parameter in system 
designe 

(2) Stochastic Inputs - Newton's method for handling 
saturation with stochastic inputs jis quite 
impressive. The method involves the derivation of 
a modified Wiener-Hopf equation with the insertion 
of a LaGrange multiplier. Since explicit solution 
of the multiplier is difficult, graphical procedures 
are usually resorted to, for determination of the 
multiplier. The method seems to have no restrictions 
as long as the use of the rms value of the 
saturating signal is an acceptable criterion for 
constraining saturation. If this is unacceptable, 
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then there seems to be no choice, analytically 
speaking, but to separate the design problem into 
two parts: a linear and a non-linear mathematical 
model, employing non-linear techniques after 
saturation occurs. 
(c) Bandwidth Minimization 
The above paragraph indicated one reason why design by 
this method called for increased bandwidth. There is 
another reason: Since the input signals are represented by 
their power-density spectra, or energy-density spectra 
(for transient inputs), minimization of ISE or MSE will 
demand that the system transmit all those frequencies in 
the input with negligible amplitude.and phase change. If 
the input signal is aperiodic, its spectrum includes all 
frequencies; if it is periodic, it includes discrete 
frequencies which are integer multiples of the fundamental 
out towards infinity. Hence, for minimum error, namely 
zero, an infinite bandwidth is called for. 
Excessive bandwidth has several undesirable features: 
(1) It allows the transmission of excessive noise; 
(2) It causes unwanted saturation by raising the 
internal signal level; along with this it 
increases the requirements for components 
operating at high power levels, since output peak 
power is generally associated with the higher 
frequencies; 
(3) It complicates the compensation, since the high 


frequency dynamics of the components cannot be 
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neglected in determining their transfer function, 
and the compensation must effectively offset these 
additional components, 

This reference presents an analytical method for 
minimizing bandwidth by transposing the problem to one of 
minimizing the mean-square value of a transmitted signal. 
The method is good for any configuration and for both tran- 
sient and stochastic signals. It uses a bandwidth test. 
For the test, a stationary stochastic noise signal 
(characterized by that to which the system is expected to be 
subjected) is used as the input to the system. The result- 
ant output is passed through a filter whose rms output is 
measured. A standard system is defined with variable 
bandwidth. It is fed the same noise signal and its output 
passed through a filter. The rms value of the control system 
filtered output is correlated with the control system 
bandwidth by comparison with the standard system filtered 
output. The bandwidth of the standard system is adjusted 
until the two rms filtered outputs are equal; by definition, 
the two bandwidths are then @qual. Only one restriction 
exists: The standard system, along with the noise source 
and filter, must produce a monotonically wieTeeoihe rms 
output with increasing bandwidths; otherwise, minimizing the 
control system filtered output would not necessarily 
minimize its bandwidth. By this scheme the problem of 
minimizing bandwidth is equivalent to determining that 


overall system transfer (or weighting) function that 
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minimizes the mean square value of its filtered output 
subject to the constraints required by the performance 
index of the normal input signal. 

It follows that additional specifications are required 
for establishing the noise source, the standard system and 
the filter. A frequently chosen noise source is white 
noise because it has components at all frequencies and is 
easily handled analytically. The filter is chosen with 
regard to the cutoff characteristics desired in the system; 
the filter must have a finite mean-squared output since this 
is the signal to be minimized. This reference recommends 
and uses a pure differentiator of (n=1)th order, requiring 
system cutoff at the rate of 1/we The standard systems 
recommended and used are the simple binomial filter or the 
Butterworth filter (see p.220, reference 2) chosen so that 
there is produced a finite value of mean-square filtered 
output. Since the problem is now transposed to that of 
minimizing the noise transmitted through the system, subject 
to the constraints imposed by the performance index, another 
problem in variational calculus results. The solution 
proceeds as indicated in the paragraph about saturation: 
The LaGrangian multiplier is introduced which eventually 
leads to a Wiener-Hopf equation. 

(d) Stability 

All of the above techniques have system stability 
inherent in them because of the "spectrum factorization" 


technique required to solve the Wiener-Hopf equation. 
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There is another matter to consider: These techniques use 
equivalent cascade compensation, requiring manipulation back 
to the feedback loop after Ga (s) has been determined. The 
possibility exists that, if the plant elements are not 
minimum phase (zeros in RHP), G_.(s) introduces a pole to 
cancel the RHP zero, and thus G.e(s) is unstable in its own 
right. To prevent this the above techniques require that 
the plant elements be stable. This is not a restriction, 
since all that is necessary is to provide some feedback 
loop around the plant elements to make them stable, and use 
the equivalent transfer function as the fixed elements for 
the rest of the design. 

(e) Disturbances 

The above line of reasoning applies to disturbances 

also: The fact that they are applied at a place other than 
the input merely requires that they be manipulated out to 
the input. This modified input is then used in the design. 
6e5 Graham & Lathrop Me thods1l 

These methods are not Analytical Design Methods, but 
their objective is the same: The Synthesis of "Optimum" 
Transient Response. This reference, first of all, seems to 
have about the most complete exposé in print on the merits 
of the various performance indices. Because of the exposé 
it is concluded that the best index is Integral-of-Time- 
Multiplied-Absolute-Error (ITAE). Since this function is 
not analytic, strictly analytical methods cannet be applied. 

An almost exhaustive study was made on the following 


criteria, as applied to a second order system with a step 


we 





input (this was a zero-displacement error system whose trans- 


fer function was of the form Cc (s) = yee) 
R S742 °S+! 


R= fettdt , B = { le)| at 
RB = (le®|de R= fe e“(t) dt 
B = ( edt Re = [tered 
PR = (+ elidt F {ew [at 


Plots of the various criteria versus damping ratio resulted 


if 


Ul 


in the following: 

(1) P, and Py dictate a 'f=O for minimum error and were 
rejected on this ground. 

(2) Pa called for ~=,5 but had very little selectivity; 
but it can be handled analytically and can be 
mechanized on a computer. 

(3) Po and Ps called for 2-53 both are easily 
mechanized for the computer; both have better 
selectivity than Pz with Ps giving much sharper 
selectivity. 

(4) Pe, PyosPi1,5 called for fx .6-.8 and gave good 
selectivity. All were rejected because they are 
too difficult to handle either analytically or by 


computer. 
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PosPa5Ps were retained and applied to a zero velocity 


error system of form 


and thence 


Cis = 2 St | 
i ) S*+ 27S! 


applied to higher order systems. It was found that P.(ITAE) 
was the only criterion which retained any selectivity. It 
was, therefore, concluded that the ITAE criterion was best. 
(It is interesting to note that one of the advantages of ISE 
criteria offered by Newton is the fact that it is not very 
selective, thus allowing wide latitude practically in 
adjusting the system.) 

Next in the study, a set of limitations is established: 
The servo must be a "duplicator"; namely, it will be 
subjected only to a step input and it will have zero-displace- 
ment error (meaning the closed loop transfer function has 
unity in the numerator). With these restrictions the servo 
is merely a low pass filter. The servo could, therefore, 
be represented by standard forms of either a binomial filter, 
or a Butterworth filter, or a "minimum ITAE" filter. It is 
a simple algebraic matter then, to force the coefficients of 
the servo characteristic equation to equal the coefficients 
of the standard forms. It was found that the minimum ITAE 
filter response combined the good response time of the 
Butterworth filter with the smaller oscillation (and over- 
shoot) of the binomial filter. However, in the attempt to 
extend the method to zero-velocity and zero-acceleration-= 
error systems, step-function responses to the minimum ITAE 


filter were no better than, and in some cases, worse than, 
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the binomial filter. 

The use of standard forms does not involve solving for 
the roots of equations, idee or graphical constructions, 
integration, or inverse Laplace transformations. As the 
reference suggests, it is a true synthesis method, in 
that it dada directly to a description of the required 
system in terms of its design parameters. There are a 
large number of servos which fall within the imposed restric- 


tions. 


6.4 Methods of Zaborszky and Diesel°*°°° 


Looking behind all the indirectness of feedback control 
system design, it develops that actual system design is based 
on the one ultimate measure of performance: How much error 
is there at the times when the system output is utilized? 
Zaborszky and Diesel attempt a mathematical formulation of 
such a measure of performance. In contrast to most standard 
techniques which concentrate on isolated phases of performance 
such as transient or steady state, this measure unites all of 
these, favoring none, in a single concepte This measure not 
only concerns itself with the amount of the error, but also 
the times when the output is being utilized. 

This measure is formulated by the following line of 
reasoning’ 

(1) The specific environment and function of each control 
system will set a specific penalty valuation on errors of a 
given size. This penalty valuation can be expressed in the 
form of a penalty function F(e), a single-valued function of 
the error e. Error itself is a function of time, and the 
valuation may vary with extraneous parameters, or with time. 
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Therefore let F(e) = F(e), oe Us, Je °° ™ UR) 


(2) A second element of performance evaluation is the 
time when the output is utilized. For simple systems,‘ all 
times of utilizing the output are generally equally probable, 
but this is generally not so in advanced or complicated 
systems. This means that the times elapsing from activation 
of the system to all times of utilization of its output can 
be arranged into a probability distribution p(t). 


The measure can be mathematically expressed as? 


5 = (Fle) ple) de 


where the symbol is the end sigma, the form of letter sigma 
used at the end of Greek words. For a deterministic input, 
the above gives the average value of the penalized error at 
all times when the output is utilized. If the environment 


consists of several inputs, or is stochastic, then 


bo 
>) = | Fee) dt 
where the bar denotes the averaging over the ensemble. If 


p(t) is independent of the input (and it frequently is), 


2) = {© pljdt 


Conceptually then, the above integrals are the average value 
of the penalized error at such times when the output is 
utilized, and are referred to as the "probabilistic error" 
or "end sigma error". As a performance index, the integral 


is called the "end sigma criterion", This measure unites in 
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a single concept the transient and steady states of operation 
as well as any intermediate states. None of these states is 
discriminated for or against, because p(t) determines the 
weight allotted to each state. The most fascinating thing 
about the above integral is its camplete generality: All of 
the other performance indices mentioned in this paper are 
simply special cases of this integral. 

It is possible to evaluate this integral in the S 
(frequency) domain if the functions involved are Laplace 
transformable. This means it can be used directly with a 
root locus plot without the necessity of going to the time 
domain. It can also be evaluated in the tim domain in the 
same manner as indicated in Chapter 4. Use of this criterion 
for optimizing a system leads again to the Wiener-Hopf 
integral equation. The second reference indicated gives a 
complete design process to obtain the impulse response function 
for the system which has the smallest average square error 
for such times when its output is used. All of the necessary 
equations are presented in an organized (albeit complicated) 
form for the programming of a digital computer. As they 
indicate, the solution of problems of a complexity common in 
control engineering practice today, requires computers. 

6.5 Methods of Schultz and Rideout’® 


These writers propose a general integral of the form 


c= (+ (ew, t] dt 


which is the same as that of Zaborszky and Diesel 
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with no consideration to the times of utilization of the 
output. However, their attitude is the recognition of the 
fact that physical systems can never respond instantaneously, 
and therefore the immediate initial error should be ignored 


by a suitable delay. They express error as follows: 


e(t,t) = nr (¢-t) — c) 


where r(t-%) represents the input delayed by an amount <r 
and c(t) is the output. This error is called "delay-error" 


and it is the one used in the indices as follows: 


bo 

E(t) = { Fle), ¢,t]dt 
Their study shows that the integral can be used with 
transient or stochastic inputs, and that an optimum value 
of * exists for a given set of parameters. They have also 
shown that IADE (Integral of Absolute value of Delay-Error), 
ISDE (Integral of Square of Delay-Error), and ITADE (Integral 
of Time-weighted Absolute value of Delay-Error) can all be 


easily instrumented on an analog computer. 
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CHAPTER 7 
CONCLUSIONS 


7ol General 

It is impossible to formulate in detail a universal 
approach to all servo design problems, but it is possible 
to list in proper sequence a number of design steps which 
can serve as a guide. As mentioned in the beginning, 
feedback control systems have penetrated all walks of life. 
All kinds of engineers and mathematicians are in the field, 
all developing methods to fit their own particular problems. 
There is another large group of persons in this field-- 
this group never really designs servos per se--they look 
to the field for stimulation and challenge to their 
ingenuity: the mathematics and graphical and block diagram 
manipulations are "fun", and a major portion of them are not 
hard to learne And so long as the de sign is a paper design, 
no "hardware" experience is really necessary. As automation 
make s further strides, less experience is needed, in one 
sense, because more and better pieces of hardware become 
available. Every day it is more nearly possible to realize 
a compensation which was not realizable a few years agoe 
But at the same time, two other things are happening? 1) 
cost control becomes more important, and 2) advanced 
systems used with refined components require a more 
sophisticated analysis. 


Given a set of specifications, many solutions to the 


he 





design problem are possible. Owing to the great variety 
of techniques and problems, experience, in the final 
analysis, plays just as important a role as the use of 
methods. As Gillet® puts it, the principle most often 
overlooked is that a feedback control system constitutes 
an entity, and each component of it must be considered as a 
part of the overall system. We would presume that this 
principle never really strikes home until one has attended 
the school of “hard knocks". 
702 Analytical Design Methods 

There can be no doubt that these techniques are rapidly 
coming into vogue. The general reason for this is that 
higher performance is continually being demanded: Feedback 
Control Systems accomplish all kinds of sophisticated jobs 
that were not done ten years ago or even conceived twenty 
years ago. The rapid development of improved components 
continues to make it easier for the designer to translate 
a complicated paper design into a useful physical system. 

The first hurdle in the use of these techniques is 
to bring oneself to accept the idea that a performance 
index can properly express the specifications of the system. 
Once this has been accepted, application of the me thods 
becomes much more palatable. The only two criteria in 
extended use are Mean Square Error (MSE) or Integral Square 
Error (ISE) and Integral Time Weighted Absolute Error 
(ITAE). Far purely analytical techniques only MSE can be 


usede 


100 


Since the objective of the design is the minimization 
of the performance index, a particular advantage of these 
methods lies in their ability to recognize inconsistent 
specifications, in the sense that if the desired performance 
index is less than the theoretical minimum value, then the 
specification cannot be fulfilled. 

Another advantage lies in the fact that the 
procedures are readily susceptible to modification to 
restrict the bandwidth or to limit saturation tendencies. 

Another advantage lies in the methods’ ability to 
handle virtually any kind of input or output. Whereas 
most Trial and Error procedures hinge on the desired output 
being equal to the input, these methods are generally 
indifferent as to what the desired output is. This means 
they are readily able to accommodate noisy signals. For 
that matter, they can handle any stationary stochastic 
signal as long as it is possible to obtain the correlation 
function, or correlation function transform. 

Another advantage, albeit a philosophical one, is that 
the mathematics associated with the methods is the same 
type of mathematics used in information and statistical 
theory and in all advanced communication theory. Persons 
schooled and trained to think in those lines can readily 
adapt their thinking to the Analytical Design Methods. 

Besides the indicated disadvantage of only one 
performance index, another disadvantage, that could be 
important sometimes, is the large number of numerical 


calculations associated with these methods. In one sense, 
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this is not serious since any method requires numerous 
calculations if the system is an advanced one; besides, 
computers are readily available today to do much of this 
detail 

From the viewpoint of minimizing MSE, these methods 
bring home very clearly three theoretical performance 
limitations on linear systems: 

1) Noise and disturbances make it impossible for a 
system to establish equality between the desired and actual 
output--which means the minimum MSE is not zero under these 
conditions, nor as low as it could be in the absence of 
the noise or disturbance. 

2) The best compensation cannot overcome the effect 
of a pure time delay in the fixed elements. This means 
that a feedback control system cannot predict the future 
value of a signal with zero error, since it inherently 
must operate on present or past information. 

3) There is no way to eliminate completely the effect 
of a non-minimum-phase fixed element. Intuitively the 
cancellation of a zero in the RHP can be approximated but 
not completely accomplished because of the threat of 
instability. Thus zero MSE cannot be obtained. 

We know that a practical performance limitation is 
that of saturation. Jn the Trial and Error Methods, 
there is a tendency to overlook this very important point 
until one attains a great deal of experience. In the 
Analytical Design Methods, saturation can be made to appear 


in the forefront from the beginning of the design. 
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705 Trial and Error Design Methods 
Cate gorically the transfer function approach will most 


directly predict frequency response while the root locus 
approach will most directly predict transient response. 

Use of the system equation most directly predicts stability. 
However, all approaches are used in actuality. 

One of the drawbacks of frequency response techniques 
is the difficulty of trying to visualize the transient 
response. The root locus method fills this gap because the 
motion of the closed-loop poles can be easily observed as 
the gain factor is varied. The root locus technique then 
serves as a fine educational tool also. Another advantage 
lies in the fact that the compensation employed or being 
investigated is readily evident and the problem of 
unrealizable compensation seldom presents itself. Direct 
determination of stability is also easy whereas the 
frequency response techniques can sometimes be misleading 
for indication of stability. 

A peculiar advantage of frequency response techniques 
(or steady state analysis) is that knowledge of the 
mathematical model is not a requirement and it is 
particularly easy to select a compensator to cause the 
response curve to take on the right shape. In root locus 
techniques (and for that matter, in Analytical Design 
Methods), the mathematical models of different parts of 
the servo must be known; in practice it may be difficult 
to determine the constants required for such a representa- 


tione 
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The major drawback to the root locus is that accurate 
plotting of the locus is a time consuming task. As 
indicated earlier, a number of theorems for approximating 
the locus are available; a number of computing schemes have 
also been developed including one, by one of the authors, 
for use on the NCR 102 computer. 

Translation from frequency response to transient 
response is also difficult at best. Here again, digital 
computers can play an important role in speeding up the 
arithmetic required in the various approximation methods. 

Analog computers of course greatly facilitate 
analysis and design. In addition, good analog computers 
are extremely valuable as an aid to remaining in the linear 
zone, or observing the effects caused by moving out of the 
linear zone. To recognize their use in this respect is 
merely to recognize that nature just is not linear. 
Probably their strongest contribution is their application 
to solution of optimization problems, after the configuration 
is fixed. 

All of the various aids and charts developed in these 
techniques come down to one points they are a scheme for 
relocating the roots of the closed-loop equation by juge 
gling the roots of the open-loop equation. In one case 
the given specifications define a closed-loop equation and 
the problem is to find an open-loop equation to fit it. In 
the second more advanced case the problem is to find an 


open-loop equation in which the plant elements fit also. 
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74 Epilogue 


The servo designer has to answer three major questions; 

a) What is the Frequency Response? 

b) What is the Transient Ragone 

c) Is it stable? | 
No single method anawereuneite questions. Careful examina- 
tion of a proponent of a particular method reveals that the 
method works well for some particular type problem, but has 
limitations when extended to some other probleme 

No matter what method is considered, the first step is 

the formulation of the problem by gathering the appropriate 
specifications and plant element descriptions. Perhaps the 
next best step is the application of Analytical Design 
Theory to the development of an appropriate formula for 
compensation. It is then necessary to make some reasonable 
approximations to reduce the complexity of this formula or 
the computations associated with it, after which the 
compensation is still unduly complex and generally unre-=- 
lizable. At this point the Trial and Error Methods should 
be injected into the amalgam to ameliorate the situation. 
Simpler forms of compensation found by these methods will 
almost always yield performance close to that determined 
theoretically. How often this is true depends upon the 
broadness of the minimum of the chosen performance index. 
Furthermore several different forms of compensation can be 
found and compared to the theoretical one for their relative 


efficiencies. The reason that this comparison is possible 
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is that it is relatively easy to optimize the parameters 
after the selection of a compensator by Analytical Design 
Methods, since the system becomes one of a fixed configura- 
tion. This can also be done with an analog computer. 
However, the practical value of optimization can be over- 
emphasized. When there are more than about two free 
parameters, optimization is "but a modern, more systematic 
variation of the very old engineering practice of 


compromise” (Quote from Gi11et®) 


Very powerful tools are available today for synthesizing 
linear feedback control systems. All are of value, providing 
the designer is well trained in their use. Apparently no 
optimum set of methods exists, although there exists a 
broad enough family of methods to arrive at "almost the best 
possible" servo for a given set of conditions. It would seem 


that further refinements of the basic methods would have to 


await additional improvement in servo components. 
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APPENDIX A 
SELECTED DEFINITIONS 


l. Dynamic Response - the output response to an input that 
is a varying function of time. 
Zo vteady-State-Response = the output response to an input 
that is constant with time. 
Se lransient Response - the time variation of one or more of 
the system outputs following a sudden change in one or more 
of the system inputs or the derivatives or integrals of the 
system inputs. A given transient response must be referred 
to the type of input that caused it. 
4. Frequency Response - the variation of the output to an 
input which is a constant-amplitude variable-frequency 
sinusoide 
oo Forced Response = the time response of an output of the 
system to an arbitrary, but completely defined, variation of 
one of the system inputs. Forced response is distinguished 
from transient response in that the input variation associated 
with the forced response of a system is considered as a 
continuous time function with no discontinuities in any of 
its derivatives. (A sinusoidal input is a special case of 
a forcing input which is isolated for special attention 
because of its theoretical importance.) 
6. Transient Test Inputs - 

ae Impulse = A unit impulse is a time function that is 


infinite at t=a and zero everywhere else. It is 
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defined as follows, where §,(t=a) is a unit impulse 


function occurring at t=a: 


ff (t-a) dt =| 


§ Silke) Fie)dt =f(a) 
Si(t-a)=o, A<tKa 


b. Step = the unit step function &., (t-a) is merely 
the integral of the unit impulse Saultoale It is 
defined as follows, where §_;(t-a) is a unit step 


occurring at t=a: - 
S.,(t-a) =| S.(K-a) dx 
— KO 


ja O, Els @ 
Lait aC 
co Ramp - The unit ramp function S,, (ta) is the 


integral of the unit step 5. (ba). The unit ramp 
is defined as follows, where SD, (ta) is a unit ramp 
occurring at t=a: 


é 
S (4-0) ={ Sulx-a) dx 


= 742 tia 
t, C7A 
d. Parabolic - The unit parabolic function S.3 (ta) 
is the integral of the unit ramp S.,(t-2)- The 
unit parabolic is defined as follows, where on 
(tea) is a unit parabolic occurring at t=a: 
5 
>. -a) = J oa lX-aldx 
©, < Gia 
Str t 7a | 
e. Displaced cosine = thé dfsplaced cosine is one 360° 
segment of a cosine wave displaced so that it is 


always positive going: 
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Te isgloce = Sle -S.,(t-a4 22})1- ae w(t-<) 


ine, 


Note that the first derivative of the displaced 
27° 


cosine is zero for times t-a and t-a TO 
Note that all of the above functions are equal to 
zero for all t4a, and that they are discontinuous, 
or one or more of their derivatives are discontinu- 
ous at the instant of occurrence. 
7. lhe Convolution Integral 

a. If y(t) is the input, x(t) the output, and w(t) the 

impulse response of the system, then the output x 


can be found by evaluating the convolution integral 


r(+) ={ w(t) u (+-t) dt, 


v(t) ={“w(t-t) 4 (4) dt 


be If the system being studied is a physical system, 


wt) =o for t<o 


and the ee ae integral reduces to: 


_ ut fu (41) 4(t-t,) dt, 
(4) = (Fw uv (t- ti) 4 (ti) dt. 


c. If y(t) and w(t) are both zero for t40, then the 


then 


convolution integral reduces to 
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ut) = (a (t,) y (t-4,) dt, 


“ule = [u(t-4) 4h) A 


8. The Fourier Transform = The Fourier transform of a 


function and its inverse are defined as follows: 


FE (4)| = F(s) = fe #(t) dt 
aN [F(s)] f(t) = ay [. "est EUs)ds 


where S = the complex variable m+ LW 

The Fourier transform is applicable to functions that exist 
for all time t. To insure the existence of the Fourier 
transform of a function, Dirichlet's conditions must be 
satisfied. 

9. The Laplace Transform - The Laplace transform of a 


function and its inverse Hs aga as follows: 


Slt] 4 Fs) & Le" ele dt 
S Tris] £44) 8 Bf e*ru)ds 


where s = the complex variableg+,W « Note that the Laplace 





transform is used for functions that are zero for t€0O. The 
constant c is used in the inverse as a convergence factor 
that enables one to apply the Laplace transform to functions 
whose Fourier transforms do not existe A function must 

also satisfy the Dirichlet conditions to be Laplace 
transformable o 

10. Gain = Gain of a system or element is the ratio of 


magnitude of the output with respect to the magnitude of 
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sinusoidal input. The frequency and conditions of operation 
and measurement must be specified. 

lle Nyquist Diagram - The Nyquist Diagram is a closed polar 
plot of a loop transfer function from which stability may be 
determined. For a single-loop system, it is a map on the 
F(s) plane of an s-plane contour which encloses the entire 
right half of the s-plane, excluding poles for the loop 
transfer function which lie on the imaginary axis. 

12. Response time - Response time is the time required for 
the output first to reach a specified value after the 
application of a step input or disturbance. 

13. Rise time = Rise time is the time required for the oute 
put to increase from one specified percentage of the final 
value to another, following the application of a step input. 
Usually the specified percentages are 10% and 90%. 

14. Settling Time - The settling time of a system or 
element is the time required for the absolute value of the 
difference between the output and its final value to become 
and remain less than a specified amount, following the 
application of a step input or disturbance. The specified 
amount is often expressed in terms of per cent of the final 
Value o 

15. Gain Margin - Gain margin is the amount by which the 
magnitude of the loop ratia of a stable system is different 
from unity at phase crossover; it is usually expressed in 
decibels. 


16. Phase Crossover - Phase crossover is a point on the 
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plot of loop ratio at which its phase angle is 180°. 

17. Phase Margin = Phase margin is the angle by which the 
phase of the loop ratio of a stable system differs from 
180°. 

18. Gain Crossover - Gain crossover is a point in the plot 
of loop ratio at which the magnitude of the loop ratio is 
unity. 

19. Loop Ratio - Loop Ratio is the frequency response of 
the primary feedback to the actuating signal. Under linear 
conditions, the ratio is expressed as GH where G represents 


the forward elements and H the feedback elements. 


112 





: 


APPENDIX B 
THE FORMULATION OF THE TRANSLATION 
AND CORRELATION FUNCTIONS 


Bel Probability Density Functions 


The analysis of stochastic signals requires the use of 
probability density functions and other statistical 
characterizations such as the average value, the root-mean- 
square value or mean-square value, and the correlation 
functions. 

The probability density functions are direct measures of 
the chance of occurrence of certain events in a process. 

The. first probability density function of a stochastic ee 
pendom) variable v(t) is denoted by 
Pi (vy 5t,) & probability that the variable has a 


value Vi at time ty 


The second probability density function is denoted by 
Po(v,,t13vote)S probability that the variable has 

a value Vv, at time t, and a value vo at time to simultaneously. 

For a stationary stochastic process (one whose statistics 


are independent of time), P,(vjt,) is independent of time ty; 


Po(v15t13Vo,te) is a function only of the time difference 
(to-t ;). (Note then, that a process can be defined with one 
less variable if it is stationary). | 
In general, the average or mean value of a stochastic 
variable v(t) is given by , 
YG) & jy P(v,t) dv 
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“he méanesquare value Me r bX 
is given by: y Xe) : ( we P(v,-t) d V 
NS 


The root-mean-square (rms) value is given by the square 
root of the mean-square value. The variance of a stochastic 


process is given by 
=i 
v- 71" 


The standard deviation @& is the square root of the variance. 


It can be expressed as follows: 
a =m, 2 Vo 
am lve-V 


Since the statistics of a stationary stochastic variable 


are independent of time, the mean value is: 


v(t) & bm +f vaat 


Tyo AT 


and the mean-square value is given by 


4G). > Fees fim Lt V7Ct) dt 


Two commonly used probability density functions are the 


normal distribution and the Poisson distributione The normal 


P(v) dv = a et Ce X) 4 


where P(v)dv is the probability of finding v between v and 


distribution is given by: 





ve dv, 


The Poisson distribution is given by 
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SNe 


P(n,at) = at) € 


n! 


where P( WM, At) is the probability of finding M events in 
a time interval At, and W is the average frequency of 


occurrence of the events. 


B-2 Correlation Functions for Stationary Stochastic Signals 





The autocorrelation function Oyy Ge of a stationary 
stochastic process v(t) is defined as the mean value of the 
product of function v at time t by the function v at time 
(t + ): 


Tt 
| Ly, 
, ) = V(t) v4) = i v(t) v(t+%) dt 
-T 

A function analogous to the autocorrelation function 
for a single signal is the cross-correlation function for a 
pair of signals. The cross correlation function Oyu (+) 
between two stationary stochastic processes v(t) and u(t) 
is defined as the mean value of the product of the function 
v at time t by the function u at time (t 4%): 

t) = Veyateea) = Bm 1 (ye ulertdt 
i = Vite CE +%) = eae v(t) u(t+%) 

—T 
Some of the useful properties of correlation functions 


ares 


(1) Dy &) i yy (-) (even function) 


(2) \Wve@)| 6 Py (o) 


™ er A, &) =~ 





e eel 


(For a stationary signal this means the autocorrela- 
tion function approaches the square of the mean 


value of the signal as 2 approaches infinity.) 


(4) Diu @) = Yay Se (not an even function) 
(5) \u (%) | ¢ VQy (0) yu (0) 
(6) bm 


aga Yuu &) =V°U 


We observe from the definition of the autocorrelation 
function, that the mean square value of the signal equals 
the value of the corresponding autocorrelation function with 
zero argument: Dy (0) = y2 
B-3 Examples of Correlation Functions 

(a) Example 1 


A common type stochastic variable used is the case 
of v(t) as a rectangular wave with values + B and -p 
and with zero crossings located at event points that 
are Poisson-distributed in time. This autocorrelation 


function can be derived as: 
Za ~ 2 |T| 
RQ, (%) = * = 


(b) Example 2 
v(t) is a rectangular wave with amplitude values 


distributed in any fashion and with zero crossings 


Pr] 


Poisson-distributed in time: 


Uy (T) = ce 


where © is the standard deviation of the amplitude 


a 
+ v(t) 


distribution and Y is the mean value of the amplitude 
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(c) Example 3 
v(t) is a train of identical finite pulses whose 


starting points are Poisson-distributed in time (known 
as "shot noise"): 


hy (t) _y(#t £(++%) At + Ta 


nO 


where f(t) is the waveform of a single pulse and 


>a) 
vey) (+) dt. (This derivation is an extension of 
hon DO 
Campbell's Theorem- see p.l02, reference 2) 
(d) Example 4 


v(t) is pure or white noise: 


en (t) = 8 é, (%) 


where & is a constant that depends on how the process 
is generated, and $e) is an impulse. For example, 
if white noise is considered as a limiting case of shot 
noise generated by exponential pulses of amplitude A, 
time constant T, and area s under the pulse, then 


Ys i ar Po 
Py aS T—> Oo 


S$ = con stant 


v= 





B-4 Correlation Function Transforms 

. Because correlation functions are completely 
defined as functions of a time variable  , they are 
Fourier transformable. By convention, = times the 


Fourier transform of a correlation function is called a 


power spectrum or a power-density spectrum. Since the 
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correlation functions involve convolution of one or more 
functions, this means that we may multiply the functions 
in the frequency plane and this is the significant advantage 


of taking their transforms. 


The power density spectrum $ (5 of a stochastic process 
is defined as 


Se 
@., (S) = = les. Ky (T) at 


The cross-power-density spectrum G,,(s) between two 


stochastic processes v(t) and u(t) is defined as 


$,,(s) 4 a ep (alse 


The inverse transformations are 


oo _ 
QyA=t fs) e7ds 
= 

| pbs st 
, = Ed 
P (T) q J Buls Ss 


Since the above integrations are along the imaginary axis 
of the s-plane, the inverse transform may be rewritten in 


terms of the real frequency W ; 
4 wt 


BO 
Py = { &, (jw) é Ww 


Recalling that the mean-square value of the signal is the 
value of the autocorrelation function with T =O » then 


oY, @s (5, (ju) dw 


Because the correlation function is an even function =e oe (a 9 
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bo 
Bx (,w) can be written as ag (Ww) 9» or YY, (2) = { Buda 
This means that the area underneath the frequency function 
B (w) over the infinite frequency range (commonly 
accepted definition of power-density spectrum) is equal to 
the mean-square value of the signal. 


Useful properties of the power spectra are 


$ vy (s) = dw (-s) (even function) 
Sy (s)"= Puy (-s) 
B-5 Translation Functions 
Since transient type signals exhibit no statistical 
properties, they are not subject to description by correla- 
tion functions. Newton introduces the use of translation 
functions (see page 51, reference 2) for this type signal. 
If x,(t) is an arbitrary transient signal, then the 
autotranslation function a, ex) is defined as: 


do 
I, @) 4 { xl) X(t+t) dt 


If x,(t) ans xo(t) are two arbitrary transient signals, 


then their cross-translation function LD yatt) is defined as; 


00 
I, (t) @ \ CX (+) Xo (44%) dt 


Although these functions characterize the signal, this 
characterization is not unique; there are a number of 
different functions that can give rise to the same trans- 
lation function. Useful properties of the translation 


function are 


+, () = =) (-*) (even function) 
Z,.(t) = 22,C%) 
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Note from the definition of the autotranslation function, 
that if T =O, ne ' ni 
7m 
A _| x7(é)d 
aL ae a (o) er {x 
where Ij represents the integral-square value of the time 
function x(t) 
B-6 Translation Function Transforms 
Because of Parseval'’s Theorem (see p.44, reference 2), 
the integral Ij, stated above, can be expressed in terms of 
its transform as 
a 
T = ae (- s) x, (s) lds 
aT 4 oe 
The Fourier transform of a translation function T,, (%) is 


defined as follows? 
PO te 
Z\2 Cs) 4 ( & T,,(t) dt 


The inverse transform is defined as 


ae 
T, (7) = amd —- S 


aeetying this definition to I, 
nN sf (FiGlds 
a = I, (0) a a 


It can be shown (see pe 57, reference 2) that 


Li (s) = x, (- $) X, CS). 


The transforms of translation functions are sometimes 





called energy density spectra as contrasted to power density 


spectra in correlation functions. It is to be noted that the 
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translation functions for many commonly encountered transient 
signals are infinite since their defining integrals do not 
convergee Particular examples are the autotranslation 
functions for the step function and the ramp. Fourier 
transforms for these functions also do not exist. In such 
situations the introduction of a convergence factor will 
frequently permit a solution to be obtained. However, in 
general, the determination of the integral square value of 

a function is usually done in the time domain unless the 


Fourier transform is a rational function. 
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